arXiv: 1504.06948vl [astro-ph.GA] 27 Apr 2015 


Astronomy & Astrophysics manuscript no. aaliaph © ESO 2015 

April 28, 2015 


Observed and predicted Lyar and UV lines 
for a sample of galaxies at redshifts z<3.7 

M. Contini^’^ 

’ Dipartimento di Fisica e Astronomia, University of Padova, Vicolo deH’Osservatorio 2.1-35133 Padova, Italy 
^ School of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 


Received 


ABSTRACT 

We explore the origin of the observed Lya and other UV lines from galaxies at z < 3.7 by detailed modelling of the spectra. The 
objects are chosen among those showing a) UV - optical - near-IR lines, b) only UV lines and c) those showing Lya in the UV and a 
few optical lines. We also present UV line predictions for a sample of galaxies in the 0.0686<z<0.8829 range. The sample of galaxies 
including Lya observations in their spectra does not show peculiar physical conditions of the emitting gas, nor abnormal element 
abundances. However, the high velocity (Vj > 1000 kms“‘) component of the emitting gas is accompanied by relatively low preshock 
densities (no~' 100-400 cm“^^) leading in some cases to broad forbidden lines. Some spectra are best reproduced by shock dominated 
models in which the photoionizing source is hidden or absent. Within more than 50 galaxies modelled in this work, only a few spectra 
from galaxies at z>2.5 correspond to a starburst temperature T, > 10^ K, similar to that found in galaxies showing some activity. 

Key words, radiation mechanisms: general — shock waves — ISM: abundances — galaxies: Seyfert — galaxies: starburst — 
galaxies: high redshift 


1. Introduction 

Lya is one of the strongest lines in the UV because H is the most 
abundant element and because corresponding to the 1S-2P tran¬ 
sition. Lya /il215 line in the far-UV is now reported by obser¬ 
vations of galaxies at high redshifts (e.g. Humphrey et al 2008, 
Dawson et al 2003, Stern et al 2002, Erb et al 2010, Fosbury et 
al 2003, Vernet et al 2000, Dey et al 2005, Norman et al 2002, 
Finkelstein et al 2011 etc). The line flux is absorbed by dust, 
depending on the environments of the Lya emitting nebula. The 
grains embedded into the gaseous clouds or distributed through¬ 
out the photon’s way to Earth can lead to a signihcant drop of the 
flux (Finkelstein et al 2011). The escape probability was treated 
by Atek et al (2014) who confirmed that the Lya escape frac¬ 
tion depends on the dust extinction, but the correlation does not 
follow the expected curve for a simple dust attenuation. They 
claim that a higher attenuation can be attributed also to a scatter¬ 
ing process and that the strength of Lya and the escape fraction 
appear unrelated to the galaxy metallicity. 

Supernova (SN) events throughout galaxies lead to specihc 
shocks in the regions surrounding the SN explosion (Heng & 
Sunyaev 2008) and close to outbursts in novae and symbiotic 
stars (Contini et al 2009). Heng & Sunyaev (2008) claim that 
charge transfer reactions between hydrogen atoms and protons 
in collisionless shocks of SN remnants produce broad Balmer, 
Lyman, and other hydrogen lines with increased Lya /Hj3 . 

Generally, photoionization by active galactic nuclei (AGN) 
or by starbursts (SB) is the main heating and ionizing mechanism 
of gas throughout the galaxy. When no photoionization source is 
seen, and much of the energy is provided by the conversion of 
kinetic energy of motion into heat, e.g. when fast moving matter 
collides with ambient ISM gas, ionization and thermal energy 


are released. They will be partly radiated as recombination and 
collisional line emission (Osterbrock 1989). 

It was suggested that high redshift galaxies originate from 
merging (e.g. Ryan et al 2008). Therefore strong shocks could 
appear in some regions of the product galaxy, as e.g. in the local 
galaxy NGC 3393 (Contini 2012). Yajima et al (2012) proposed 
a formation model from major mergers of gas rich galaxies emit¬ 
ting the Lya line. They claim that at z> 2 the merger rate is 
higher in dense region and the progenitors are more gas rich. 
Lyman break technique could be employed to explore galax¬ 
ies at 0.8 < z < 2.5 (Oteo et al. 2002). Merging galaxies emit 
Lya as well as many other lines. Observations of UV lines are 
now available from galaxies at high z, in particular NV 1240, 
CIV 1550, Hell 1640, etc. and eventually OVI 1034. They are 
emitted by high ionization-level ions that are easily produced 
in shock heated gas corresponding to suitable velocities. High 
ionization-level lines are strong from C, N, O, Ne, etc. in the 
IR, optical, and UV frequency ranges. Most of the UV lines 
(e.g. OIV, OV, NV) are permitted, while e.g. [NV], [OV], etc in 
the optical and IR are forbidden and thus collisionally deexcited 
at high densities.. Permitted lines originate from recombination 
and depend on radiative processes. Both radiation and collisional 
processes should be considered in the calculation of line spectra 
rich enough to constrain the models. Composite models (shock 
-Hphotoionization) were adopted to explain the spectra from local 
merging galaxies (e.g. NGC 7212, NGC 3393, NGC 6240, Arp 
220, Contini 2013a and references therein) and were used in the 
spectral line analysis of galaxies at a relatively high z (Contini 
2014a and references therein.) 

In this paper we explore the origin of the observed Lya and 
other lines in the UV from galaxies at z < 3.7 by detailed mod¬ 
elling of the spectra. Our aim is to hnd out not only the gas phys¬ 
ical conditions and element abundances throughout the emitting 
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Table 1. The symbols in Fig. 1 


Vs 

klTlS“^ 

no 

logT. 
[10^ K] 

logD 

[cm] 

symbols 

AGN too 

100 

- 

17 

black dotted 

100 

100 

- 

19 

black dash-dotted 

200 

200 

- 

18 

blue dotted 

200 

200 

- 

19 

blue dash-dot 

100 

300 

- 

17 

magenta dotted 

100 

300 

- 

19 

magenta dash-dot 

300 

300 

- 

17 

green dotted 

300 

300 

- 

19 

green dash-dot 

500 

300 

- 

18 

red dotted 

500 

300 

- 

19 

red dash-dot 

SB 100 

100 

4.0 

19 

red solid 

100 

100 

4.7 

19 

red dash-dot 

100 

100 

5.0 

19 

red dash 

200 

200 

4.0 

19 

blue solid 

200 

200 

4.7 

19 

blue dash-dot 

200 

200 

5.0 

19 

blue dash 

300 

300 

4.0 

19 

magenta solid 

300 

300 

4.7 

19 

magenta dash-dot 

300 

300 

5.0 

19 

magenta dash 

500 

300 

4.0 

19 

black solid 

500 

300 

4.7 

19 

black dash-dot 

500 

300 

5.0 

19 

black dash 


clouds, but to determine the nature of the radiation source, such 
as a AGN, a SB or only shocks. The objects were chosen among 
those showing a) UV - optical - near-IR lines (Humphrey et al, 
Dawson et al. Stern et al, Erb et al and Fosbury et al), b) only 
UV lines (Vernet et al, Dey et al and Norman et al) and c) those 
showing Lya in the UV and a few optical lines (Finkelstein et al). 
We also present UV line predictions for a sample of galaxies (Fy 
et al 2014) in the 0.0686<z<0.8829 range, where the physical 
conditions of the emitting gas are constrained by the observed 
optical lines, in particular [OIII]4363. The total sample which 
includes galaxies in an all-sky radio survey, hard X-ray emit¬ 
ting sources, a type 2 QSR identified as a X-ray source, a type 2 
AGN at high z (3.7), a very large spacially extended Fycr emit¬ 
ting nebula, objects in a multi-band imaging survey for Fyman 
break galaxies, in a slitless spectroscopic survey for FAEs, etc. is 
rather heterogeneous. However, the observed spectra from high 
z objects accounting for the Fycr line and for enough lines con¬ 
straining the models, are few, so we gathered them from different 
classes of objects. 

We calculate for each galaxy the gas density, temperature, 
shock velocity, etc., as well as the element abundances and the 
photoionizing source. The presence of an AGN or/and of a SB, 
collision heating and ionization of the ISM gas by shocks, the 
element abundances and dust survival throughout the clouds are 
important issues for middle and high redshifts galaxies which are 
tracers of star formation and evolution. The calculation models 
are briefly described in Sect. 2. Modelling results of the galaxy 
spectra are presented in Sect. 3. Results and concluding remarks 
are discussed in Sect. 4. 


2. The calculation code 

The line and continuum ^ectra emitted by gas and dust are cal¬ 
culated by the code sum4j. The code simulates the physical con¬ 
ditions in an emitting gaseous cloud under the coupled effect of 
photoionization from an external radiation source and shocks. 
The line and continuum emission from the gas are calculated 


' http://wise-obs.tau.ac.il/~marcel/suma/index.htm 


consistently with dust-reprocessed radiation in a plane-parallel 
geometry. 

2.1. Input parameters 

To calculate line and continuum fluxes emitted from a gas neb¬ 
ula, the temperature, density and element abundances should be 
known. In a shock dominated hydrodynamical regime the den¬ 
sity profile throughout the nebula is shaped by the shock, while 
radiative and collisional heating and ionization of the gas are 
due to radiation from the external source (AGN, SB) and by the 
shock, respectively. The input parameters which characterise the 
shock are roughly suggested by the data, e.g. the shock velocity 
Vs by the FWHM of line profile and the pre-shock density no by 
the characteristic line ratios and by the pre-shock magnetic field 
Bq. We adopt Bq = 10“^ G, which is suitable to the NFR of AGN 
(Beck 2011). 

The ionizing radiation from an external source is character¬ 
ized by its spectrum and by the flux intensity. The flux is calcu¬ 
lated at 440 energies, from a few eV to keV. If the photoioniza¬ 
tion source is an active nucleus, the input parameter that refers 
to the radiation field is the power-law flux from the active center 
F in number of photons cm“^ s ' eV ' at the Fyman limit with 
spectral indices a(/y=-1.5 and ax--0.7. It was found by mod¬ 
elling the spectra of many different AGNs that these indices were 
the most suitable, in general, (see, e.g. Contini & Aldrovandi 
1983, Aldrovandi & Contini 1984, Rodrfguez-Ardila et al 2005, 
and references therein). In particular, for all the models pre¬ 
sented in the following, we use auv - -1-5, recalling that the 
shocked zone also contributes to the emission-line intensities, so 
that our results are less dependent on the shape of the ionizing 
radiation. The power-law in the X-ray domain was found flatter 
by the observations of local galaxies (e.g. Crenshaw et al 2002, 
Turner et al 2001). F is combined with the ionization parameter 
U by : t/= (F/{n c (auv -1)) ((£’/r)~“''''’'^^ - (£’c)~“''''^*) (Contini 
& Aldrovandi, 1983), where £// is H ionization potential and Ec 
is the high energy cutoff, n the density, auv the spectral index, 
and c the speed of light. If stars with a colour temperature T, are 
the photoionization source, the number of ionizing photons cm“^ 
s“' produced by the hot source is By/hv dv, where vq = 

3.29 10'^ Hz and By is the Planck function. The flux from the 
star is combined with (7 and n by /V (r/R)^=Unc, where r is the 
radius of the hot source, R is the radius of the nebula (in terms 
of the distance from the stars). Therefore, T, and t/ compensate 
each other, but only in a qualitative way, because the black body 
radiation (which depends on T* ) determines the frequency dis¬ 
tribution of the primary flux, while t/ represents the number of 
photons per number of electrons reaching the nebula. The choice 
of T, and t/ is made by the fit of the line ratios. 

In addition to the radiation from the primary source, the ef¬ 
fect of the diffuse radiation created by the hot gas is also cal¬ 
culated, using 240 energies for the spectrum. The geometrical 
thickness of the emitting nebula D determines whether the model 
is radiation-bounded or matter-bounded. The dust-to-gas ratio 
(d/g) and the abundances of He, C, N, O, Ne, Mg, Si, S, A, Fe 
relative to H, are also accounted for. We adopt an initial grain 
radius agr-lpm . 

2.2. Calculation process 

The code accounts for the direction of the cloud motion rela¬ 
tive to the external photoionizing source. A parameter switches 
between inflow (the radiation flux from the source reaches the 
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Fig. 1. Left; the Lya /Hy6 ratio calculated by the grid of models for the NLR of AGN (symbols are given in Table 1); centre ; for SB 
(symbols in Table 1) right; for shock dominated spectra (red line ; D=10'^ cm; blue line ; D=10*® cm) 


Table 2. Modelling the spectra observed by Humphrey et al (2008) 



0211-122 * 

FWHM* 

MHl 

MH2 

MH3 

av 

0406-244 * 

MH4 

MH5 

av 

z 

SFR(Mo yr-*) 

Lya 1216 

2.340 

24.85 

0.21+0.01 

1000 

9.7 

13.3 

9.9 

10.0 

2.440 

911. 

7.0+0.4 

25.6 

21.2 

25.5 

NV 1240 

0.35+0.03 


1. 

3.87 

0.005 

0.2 

<0.2 

0.005 

10.9 

0.25 

SilV+OIV 1403 

0.074+0.008 


0.33 

3.33 

0.01 

0.15 

0.3+0.1 

0.04 

10.0 

0.27 

NIV]1485 

0.024+0.005 


0.07 

0.83 

0.007 

0.04 

<0.1 

0.02 

2.17 

0.069 

CIV 1549+ 

0.65+0.02 


1. 

7.67 

0.15 

0.5 

0.9+0.1 

0.31 

20.7 

0.772 

Hell 1640 

0.36+0.01 

600 

0.83 

0.67 

0.4 

0.42 

1.2+0.1 

1.44 

5.27 

1.53 

OIll]1667 

0.028+0.003 


0.13 

1.23 

0.04 

0.1 

0.2+0.1 

0.09 

3.6 

0.17 

CIll] 1909 

0.17+0.03 


0.33 

3.67 

0.46 

0.5 

0.9+0.1 

0.8 

9.9 

I.O 

CIl] 2326 

- 


0.2 

1.57 

0.1 

0.17 

0.6+0.1 

0.46 

3.85 

0.54 

[NeIV]2423 

0.11+0.01 


0.07 

0.7 

0.03 

0.06 

1.0+0.3 

0.09 

8.0 

0.3 

[NeV] 3426 

<0.3 


0.07 

0.7 

0.02 

0.05 

0.5+0.1 

0.05 

1.8 

0.1 

[01113728+ 

0.39+0.08 

730 

0.17 

1.33 

0.45 

0.48 

2.6+0.2 

2.65 

3.7 

2.67 

[NeIII]3869 

0.31+0.08 

390 

0.13 

0.67 

0.34 

0.36 

1.0$+0.1 

1.13 

2.0 

1.15 

[011114363 

<0.4 


0.02 

0.2 

0.032 

0.04 

- 

- 

- 

- 

Hell 4686 

<0.2 


0.07 

0.07 

0.03 

0.036 

0.21+0.06 

0.2 

0.4 

0.21 

quiescent 

<0.5 


0.33 

0.33 

0.33 

0.33 

I.+0.7“* 

1. 

1. 

1. 

[OIII] 5007+ quiescent 

2.81+0.1 

< 500 

0.4 

2.67 

4.3 

4.2 

12.6+0.1 

13. 

10.0 

12.9 

[Ollll 5007+ perturbed 

1.61+0.06 

1300 

- 

- 

- 

- 

- 

- 

- 

- 

[Oil 

<0.2 


0.10 

0.67 

0.09 

0.11 

- 

1.3 

0.58 

1.28 

Ho- quiescent 

0.77+0.15* 

< 500 

1 

1 

1 

1.0 

3.3 

2.9 

3.1 

2.9 

Hct perturbed 

0.23+0.08 

- 

- 

- 

- 

- 

- 

- 

- 

- 

[Nlll 6585+ 

<0.2 

- 

0.43 

1.33 

0.58 

0.6 

- 

- 

- 

- 

[Sill 6725+ 

<0.2 

- 

0.013 

0.13 

0.05 

0.055 

- 

- 

- 

- 

Vs (kms^*) 

<0.2 

- 

1000 

500 

300 

- 

- 

100 

700 

- 

no (cm“^) 

- 

- 

400 

300 

300 

- 

- 

100 

300 

- 


- 

- 

- 

- 

20 

- 

- 

1.3 

- 

- 

D® 

- 

- 

1 

1 

0.01 

- 

- 

1 

0.028 

- 

C/H® 

- 

- 

3.3 

3.3 

3.3? 

- 

- 

3.3 

3.3 

- 

n/h’ 

- 

- 

1 

1 

1 

- 

- 

1. 

1. 

- 

0/H’ 

- 

- 

6.6 

6.6 

6.6 

- 

- 

6.6 

6.6 

- 

s/h'* 

- 

- 

- 

- 

- 

- 

- 

0.02 

0.02 

- 

HjS calc* 

- 

- 

0.1 

0.037 

1.1 

- 

- 

0.05 

0.006 

- 

w 

- 

- 

0.014 

0.58 

0.40 

- 

- 

0.833 

0.167 

- 


* in 10 *®ergcm^s *;^inkms ^ flux of narrow Hn is 6.4±0.5 10 *®ergcm^s *; ^ flux of H/3 is 3.3±0.2 10 ^ in 10*® photons cm ^ s ' eV 'atthe 

Lyman limit; ® in 10*® cm; ^ in 10“^ units; ® in erg cm“^ s“* 


shock front edge of the cloud) and outflow (the flux reaches the 
edge opposite to the shock front). The calculations start at the 
shock front where the gas is compressed and thermalized adia- 
batically, reaching the maximum temperature in the immediate 
post-shock region (T~ 1.5 x 10^ (Vs /lOO km s“*)^). T decreases 
downstream by the cooling rate and the gas recombines. The 
downstream region is cut into a maximum of 300 plane-parallel 
slabs with different geometrical widths calculated automatically, 
in order to account for the temperature gradient (Contini 2009 
and references therein). In each slab, compression (n/no) is cal¬ 
culated by the Rankine-Hugoniot equations for the conservation 
of mass, momentum and energy throughout the shock front (Cox 


1972). Compression ranges between 4 (the adiabatic jump) and 
> 100, depending on Vs and Bq. The stronger the magnetic held, 
the lower the compression downstream, while a higher shock ve¬ 
locity corresponds to a higher compression. The cooling rate is 
calculated in each slab by free-free (bremsstrahlung), free-bound 
and line emission. Therefore, most emission lines must be calcu¬ 
lated in each slab even if only a few ones are observed because 
they contribute to the temperature slope downstream. 

In pure photoionization models, the density n is constant 
throughout the nebula. In models accounting for the shocks, both 
the electron temperature and density ng show a characteris¬ 
tic prohle throughout each cloud. The density reaches its upper 
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Table 3. Modelling the spectra by Humphrey et al (2008) 


flux* 

0529-549 

MH6 

MH7 

0828 

MH8 

MH9 

av 

1138-262 

MHIO 

MHll 

4C+23.56 

MH12 

MH13 

z 

SFR(Mo yr-') 
Lya 1216 

2.575 

150 

27.0 

21.2 

2.572 

86. 

21.7 

45. 

32.5 

33.8 

2.156 

56. 

29. 

27.6 

2.479 

163. 

2.17 

47.5 

30.4 

NV 1240 

- 

0.07 

10.9 

0.9 

24. 

0.02 

2.6 

- 

0.024 

0.017 

0.31 

0.03+ 

0.06+ 

SilV+OIV 1403 

- 

0.09 

10. 

0.78 

2.3 

0.05 

0.3 

- 

0.02 

0.013 

0.05 

0.19 

0.13 

NIV]1485 

- 

0.02 

2.17 

0.18 

1.3 

0.02 

0.16 

- 

0.002 

0.002 

0.06 

0.08 

0.02 

CIV 1549+ 

- 

0.2 

20.7 

4.5 

30.3 

0.2 

3.5 

- 

0.054 

0.01 

0.51 

1.63 

0.33 

Hell 1640 

- 

0.83 

5.27 

2.57 

5.8 

1. 

1.5 

- 

0.6 

0.011 

0.43 

3.7 

1.32 

Om]1667 

- 

0.065 

3.6 

0.46 

3.15 

0.1 

0.43 

- 

0.04 

0.009 

0.11 

0.37 

0.11 

CIII] 1909 

- 

0.55 

9.9 

1.4 

4.2 

1.1 

1.43 

- 

0.46 

0.1 

0.28 

0.7 

1. 

CII] 2326 

- 

0.56 

3.85 

0.48 

3.4 

0.9 

1.17 

- 

0.35 

0.1 

0.22 

0.19 

0.85 

[NeIV]2423 

- 

0.04 

8. 

1.2 

2. 

0.05 

0.26 

- 

0.02 

0.0027 

0.37 

0.26 

0.06 

[NeV] 3426 

<2.1 

0.018 

1.8 

0.76 

8.85 

0.012 

0.98 

- 

0.01 

0.0028 

0.74 

0.2 

0. 

[OII]3728+ 

3.2 

3.3 

3.7 

1.75 

2.1 

3.7 

3.5 

1.44 

1.35 

1.2 

2.57 

2.7 

3.4 

[NeIII]3869 

< 2.1 

1.0 

2. 

- 

1.6 

1.48 

1.5 

0.82 

0.68 

0.29 

0.74 

1.36 

1.5 

[OIII] 4363 

<1 

0.04. 

0.63 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Hell 4686 

<1 

0.12 

0.4 

0.29 

0.44 

0.15 

0.18 

<0.88 

0.09 

0.002 

<0.86 

0.53 

0.2 

W 

1. 

1. 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

[OIII] 5007+ 

10.7 

9.0 

10 

9.53 

8. 

11.8 

11.4 

6. 

6.3 

6.16 

13.05 

12.5 

13.5 

[OI] 6300+ 

0.58 

1.1 

- 

- 

- 

- 

- 

<0.59 

0.5 

0.002 

<0.28 

0.7 

1 

Ha 

5.2 

3.0 

3.1 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

[Nil] 6585+ 

2.6 

2.76 

0. 

- 

- 

- 

- 

2.35 

1.9 

0.8 

0.8 

1.0 

1.86 

[SII] 6725 

1.58 

0.6 

0. 

- 

- 

- 

- 

0.85 

0.63 

0.1 

0.86 

0.8 

0.7 

Vs (kms^') 

- 

180 

700 

- 

700 

1200 

- 

- 

300 

300 

- 

200 

200 

no (cm“3) 

- 

100 

300 

- 

400 

400 

- 

- 

300 

300 

- 

300 

300 

T, 2 

- 

- 

- 

- 

- 

- 

- 

- 

- 

4 

- 

25. 

- 

U 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.8 

- 

0.16 

- 


- 

1.0 

- 

- 

- 

6 

- 

- 

8 

- 

- 

- 

4 

O'* 

- 

1 

0.028 

- 

0.019 

0.3 

- 

- 

0.01 

0.01 

- 

1. 

0.01 

C/H^ 

- 

3.3 

3.3 

- 

3.3 

3.3 

- 

- 

3.3 

3.3 

- 

0.8 

0.8 

N/H 5 

- 

1 

1 

- 

1 

1 

- 

- 

1. 

1 

- 

0.5 

0.5 

O/H^ 

- 

6.6 

6.6 

- 

6.6 

6.6 

- 

- 

6.6 

6.6 

- 

6.6 

6.6 

S/H^ 

- 

- 

- 

- 

- 

- 

- 

- 

0.3 

0.3 

- 

0.3 

0.3 

np calc 

- 

0.06 

0.0058 

- 

0.0089 

0.29 

- 

- 

0.62 

1.6 

- 

0.59 

0.14 

W 

- 

- 

- 

- 

0.8 

0.2 

- 

- 

- 

- 

- 

- 

- 

' in 10" 

ergcm“^ s“ 

‘;2in 10"* 

K; 3 in 

10*® photons 

cm“^ s“ 

* eV"' 

at the Lyman limit; ^ in 

10*® cm; 

5 in lO""* 

units; ^ in erg cm ^ s 



limit at a certain distance from shock-front and remains nearly 
constant, while ng decreases following recombination. 

The primary and secondary radiation spectra, calculated by 
radiation transfer, change throughout the downstream slabs, each 
of them contributing to the optical depth. In each slab of gas 
the fractional abundance of the ions of each chemical element 
is obtained by solving the ionization equations which account 
for photoionization (by the primary and diffuse secondary radi¬ 
ations and collisional ionization) and for recombination (radia¬ 
tive, dielectronic), as well as for charge transfer effects, etc. The 
ionization equations are coupled to the energy equation when 
collision processes dominate (Cox 1972) and to the thermal bal¬ 
ance if radiative processes dominate (Williams 1967). The latter 
balances the heating of the gas due to the primary and diffuse ra¬ 
diations reaching the slab with the cooling due to line emission, 
dust collisional ionization and thermal bremsstrahlung. The line 
intensity contributions from all the slabs are integrated through¬ 
out the cloud. In particular, the absolute line fluxes referring to 
the ionization level i of element K are calculated by the term 
tiKii) which represents the density of the ion X(i). We consider 
that n/f(i)=X(i)[K/H]nH, where X(i) is the fractional abundance 
of the ion i calculated by the ionization equations, [K/H] is the 
relative abundance of the element K to H and hh is the density 
of H (by number cm“^). In models including shock, hh is calcu¬ 
lated by the compression equation in each slab downstream. So 
the element abundances relative to H appear as input parameters. 
To obtain the N/H relative abundance for each galaxy, we con¬ 
sider the charge exchange reaction N^-i-H ^ Nh-H^ (Steigman 
et al. 1971). Charge exchange reactions occur between ions with 


similar ionization potential (I(H^)=13.54 eV, I(N^)=14.49 eV 
and I(0^)=13.56 eV). It was found that N ionization equilib¬ 
rium in the ISM is strongly affected by charge exchange. This 
process as well as O^-i-H ^ O+H^ (Field & Steigman 1971) are 
included in the SUMA code. The N^/N ion fractional abundance 
follows the behaviour of 0^/0 so, comparing the [NII]/Hj6 and 
the [OII]/Hy6 line ratios with the data, the N/H relative abun¬ 
dances can be easily determined (see Contini et al. 2012). 

Dust grains are coupled to the gas across the shock front 
by the magnetic field. They are heated radiatively by photoion¬ 
ization and collisionally by the gas. up to the evaporation tem¬ 
perature (Tdust ^ 1500 K). The distribution of the grain radii 
downstream is determined by sputtering (Dwek 1981), which de¬ 
pends on the shock velocity and on the gas density. Throughout 
shock fronts and downstream, the grains might be completely 
destroyed by sputtering. The dust-to gas ratio djg affects the mu¬ 
tual heating and cooling of gas and dust. 

The calculations proceed until the gas cools down to a tem¬ 
perature below 10^ K (the model is radiation bounded) or the 
calculations are interrupted when all the lines reproduce the ob¬ 
served line ratios (the model is matter bounded). In case that 
photoionization and shocks act on opposite edges, i.e. when the 
cloud propagates outwards from the radiation source, the cal¬ 
culations require some iterations, until the results converge. In 
this case the cloud geometrical thickness plays an important role. 
Actually, if the cloud is very thin, the cool gas region may dis¬ 
appear leading to low or negligible low ionization level lines. 

Summarizing, the code starts in the first gas slab adopting 
the input parameters given by the model. Then, it calculates the 
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density, the fractional abundances of the ions from each level for 
each element, free-free, free-bound and line emission fluxes. It 
calculates Teby thermal balancing or the enthalpy equation, and 
the optical depth of the slab in order to obtain the primary and 
secondary fluxes by radiation transfer for the next slab. Finally, 
the parameters calculated in slab i are adopted as initial condi¬ 
tions for slab i-i-1. Integrating the line intensities from each slab, 
the absolute fluxes of the lines and of bremsstrahlung are ob¬ 
tained at the nebula. The line ratios to a certain line (generally 
HjS for the optical-UV spectrum) are then calculated and com¬ 
pared with the observed data, in order to avoid problems of dis¬ 
tances, absorption, etc. The number of the lines calculated by 
the code (more than 300) does not depend on the number of the 
observed lines nor on the number of input parameters, but on the 
elements composing the gas. 

2.3. Selection of the models. The grids 

The physical parameters are combined throughout the calcula¬ 
tion of forbidden and permitted lines emitted from a shocked 
nebula. The ranges of the physical conditions in the gas are de¬ 
duced, as a first guess, from the observed line ratios because they 
are more constraining than the continuum SED. We make use 
of the grids presented by Contini & Viegas (2001a,b) for AGN 
(GRIDl) and SB (GRID2), respectively for a first estimate of the 
gas physical conditions. The grids were calculated by a gradual 
increase of the input parameters adapted from the observed is¬ 
sues. Grid results give a rough approximation to the observed 
spectra. Diagnostic diagrams for Lya fli/3 calculated by the code 
SUMA for ejected clouds (Fig. 1) show that radiation dominated 
models (H-shocks) lead to Lya /H/3 between 20 and 30, with some 
higher values (up to 50) for low velocity, high density extended 
clouds. For shock dominated clouds Lya fli/3 ranges between 
-150 and-15. 

The parameters are then refined by the detailed modelling of 
the spectrum. Specific grids of models are calculated for each 
galaxy modifying the parameters on small scales in order to re¬ 
produce as close as possible all the observed line ratios in each 
spectrum. We generally consider that the observed spectrum is 
’’satisfactorily” fitted by a model when the strongest lines are 
reproduced by the calculation within 20% and the weak ones 
within 50%. ’’Satisfactorily” means that the approximation can 
be accepted 1) in the context of limits set by the author com¬ 
munity, 2) considering the ensemble of all the line fit, 3) within 
the limits of the observational error, 4) if it can be improved only 
adopting unsuitable gas conditions, and 5) because, nevertheless, 
adding some information. 

The final gap between observed and calculated line ratios is 
due to the observational errors both random and systematic, as 
well as to the uncertainties of the atomic parameters adopted by 
the code, such as recombination coefficients, collision strengths 
etc., which are continuously updated, and to the choice of the 
model itself. 

The observed spectra provide averages throughout the whole 
galaxy, particularly at high redshifts. When the observed line en¬ 
semble is not satisfactorily reproduced by one model, the contri¬ 
bution of other regions within the galaxy at different physical 
conditions must be accounted for, leading to pluri-cloud models. 
The different results are summed up by relative weights (see, e.g. 
Rodriguez-Ardila et al. 2005) when they compensate each other 
(e.g. in Table 2). When they show similar results they are left as 
alternative proposals (e.g. Table 3). 

The models selected by the fit of the line spectrum are cross¬ 
checked by fitting the continuum SED. In the UV range the 


bremsstrahlung from the nebula is blended with black body 
emission from the star population background. The maximum 
frequency of the bremsstrahlung peak in the UV - X-ray domain 
depends on the shock velocity. In the IR range dust reprocessed 
radiation is generally seen (Contini, Viegas & Prieto 2004). In 
the radio range synchrotron radiation by the Fermi mechanism 
at the shock front is easily recognized by the slope of the SED. 

The set of the input parameters which leads to the best fit 
of the observed line ratios determines the physical and chem¬ 
ical properties of the emitting gas. They are considered as the 
’’results” of modelling. 

3. Modelling single galaxy spectra 

3.1. Spectra of the radio galaxies at z- 2.5 from the 
Humphrey et al (2008) survey 

The long-slit NIR spectra obtained by the Infrared Spectrometer 
And Array Camera (ISAAC) instrument at VLT and com¬ 
pleted with optical lines (by Keck and Subaru) contain a col¬ 
lection of lines in the wavelength range 1200-7000 A. The ob¬ 
served Lya /H/3 line ratios bridge between UV and optical data. 
Humphrey et al, modelling the spectra by pure radiation and 
pure shock models, concluded that the galaxies are most prob¬ 
ably AGN in different physical conditions, with roughly solar 
abundances. We have selected dereddened (Ha /H/3 - 3) spec¬ 
tra of galaxies containing the H/3 flux in order to compare the 
Lya /H/3 line ratios (calculated at the nebula) with the observa¬ 
tion data. We wish to investigate whether an SB contributes to 
some of the observed line ratios. 

We consider all the line ratios in each spectrum (Tables 2 
and 3), trying to reproduce them by models which account for 
the coupled effect of photoionization and shocks. The models 
are described in the bottom of the tables. Models corresponding 
to F=0 are shock dominated, i.e. the photoionising flux from the 
source outside the emitting nebula is hidden or absent. Similar 
shock dominated models are used to model supernova remnant 
spectra. The FWHM of the line profiles gives a hint about the 
shock velocity. Different FWHM of the lines in the UV and in 
the optical range indicate that the spectra emitted from different 
clouds within the galaxy should be considered. So pluri-cloud 
models are adopted. The weighted sum of the line ratios calcu¬ 
lated by different theoretical models for each spectrum (e.g. col. 
7 and 11 in Table 2) is compared with the observed line ratios. 
The weights adopted to sum the models are reported in the bot¬ 
tom row of the table. 

In Table 2 we show the modelling of galaxies 0211-122 and 
0406-244. The line fluxes are referred to Ha =1 for 0211-122 
and to H/3 = 1 for 0406-244, following the Humphrey et al. no¬ 
tation. The H/3 flux is given as an upper limit for 0211-122. 
Considering that the calculated Ha /H/3 line ratio is - 3 and the 
lines are referred to Ha =1, H/3 results -0.33, which is within 
the upper limit (0.5). Moreover, the spectrum is rich in lines 
from different elements in different ionization levels which can 
constrain the models. We focus on the [OII]3727-i-/[OIII]5007-i- 
and CIV1549 h-/CIII]1909-h/CII]2326-h line ratios (The -H indi¬ 
cates that the doublet is summed) We could not find any model 
showing Lya /Ha as low as 0.2. By a shock dominated model 
with Vs =2000 kms"^and no=1300 cm“^ it is possible to ob¬ 
tain Lya fH/3 -1 (as for the 0406-244 galaxy) and even less. But 
in this case Hell 1640/HjS =9 and all the other line ratios =0 
destroy the good fit obtained by the radiation dominated model. 
High FWHM are adapted to the AGN broad line-emission region 
(BLR) (> 1000 km s ') rather than to the narrow line-emission 
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Table 4. Modelling UV spectra by Dawson et al (2003), Stern et al (2002) and Erb et al (2010) 



HDF‘ 

FWHM^ 

MDS 

CXO* 

FWHM’ 

MSS 

MSR 

8X418** 

FWHM 

ME 1-0 

ME2-0 

ME3 

MESB 

z 

2.011 



3.288 




2.3 






SFR(Mo yr-‘) 

4.73 



18.75 




17.5 






OVI 1035 

- 

- 

44.4 

1.5: 

2640 

12.8 

0.25 

- 

- 

- 

- 

- 

- 

OV 1215 

- 

- 

2.4 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Lyo- 1216 

1.68 

1270 

1.5 

18.9 

1520 

29.9 

35. 

11.27 

840 

35. 

34.5 

29.6 

26.6 

NV 1240 

2.3 

2100 

2.2 

0.6: 

1820 

0.65 

0.16 

- 

- 

- 

- 

- 

- 

SilV+OIV 1403 

- 

- 

- 

0.4 

1320 

0.3+0.86 

0.22 

- 

- 

- 

- 

- 

- 

CIV 1549+ 

2.47 

1300 

2.2 

3.5 

1350 

3.5 

3.4 

- 

- 

- 

- 

- 

- 

Hell 1640 

0.45 

1400 

1.0 

1.7 

940 

2. 

3.0 

0.31 

612 

0.3 

0.31 

0.89 

0.003 

OIII] 1667 

- 

- 

- 

0.9 

1290 

0.3 

0.4 

0.23 

235 

0.033 

0.06 

0.14 

0.02 

CIII] 1909 

0.19 

900 

0.3 

2.1 

1090 

2. 

1.9 

0.54 

225 

0.16 

0.27 

0.88 

0.2 

NelV] 2424 

0.3 

1470 

0.20 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

[Oil] 3727 

- 

- 




- 


<0.73 

-225 

0.045 

0.09 

0.65 

0.66 

US 4861 

- 

- 

0.1 

1.2 

170 

1 

I 

1 

<102 

1 

1 

1 

1 

[OIII] 5007+ 

- 

- 

0.7 

18.1 

300-430 

2.7 

17.1 

8.5 

140 

0.2 

0.49 

8. 

8.88 

[Nil] 6584+6548 

0.96 

380 

0.02 

- 

- 

- 

2.74 

<0.12 

- 

0.28 

0.48 

0.2 

0.22 

Hof n 

0.7 

240 

- 

- 

- 

- 

3. 

3.08 

155 

3. 

2.9 

3.3 

3. 

Hff b 

2.3 

2500 

0.278 

- 

- 

2.7 

- 

- 

- 

- 

- 

- 

- 

Vs (kms^*) 

- 

- 

1300 

- 

- 

1500 

200 

- 

- 

840 

600 

260 

260 

no (cm“3) 

- 

- 

300 

- 

- 

110 

350 

- 

- 

3000 

3000 

200 

200 

F^ 

- 

- 

0 

- 

- 

- 

100 

- 

- 

50 

50 

30 

- 

'j’ 6 

- 

- 

2.5 

- 


- 

- 

- 

- 

- 

- 

- 

6. 

u 

- 

- 

8 

- 


- 

- 

- 

- 

- 

- 

- 

0.7 

D (I0''’cm) 

- 

- 

9.0 

- 

- 

41.5 

100 

- 

- 

3.6 

2. 

400 

400 

C/H’ 

- 

- 

2.6 

- 

- 

3.3 

3.3 

- 

- 

3.3 

3.3 

3.3 

3.3 

N/H’ 

- 

- 

1. 

- 

- 

0.4 

0.1 

- 

- 

0.3 

1.0 

0.3 

0.3 

O/H’ 

- 

- 

6.6 

- 

- 

3.6 

6.6 

- 

- 

6.6 

6.6 

6.6 

6.6 

UJ3 calc** 

- 

- 

0.0085 

-- 

- 

0.035 

1.77 

- 

- 

10.8 

4.5 

1.65 

0.91 


' HDFx28 (Dawson et al 2003) in lO^'’ ergcm^^ g-*; Un kms-‘; ’ CXO J084837.9+445352 (Stern et al 2002) in 10^*'' ergcm^^ g-*; BX418 (Erb et al 2010) in 
2.6 10“’^ erg ® in 10*** photons s“* eV“* at the Lyman limit; ® in 10“* K ^ in 10“^ ; * in erg cm“^ s“*; 


region (NLR). Forbidden lines such as [Oil] 3727 are generally 
collisionally deexcited in the BLR due to the high density. Low 
Lya /Ha perhaps indicates strong absorption by dust in the neb¬ 
ula. The results presented in Table 2 show that the lines come 
from different clouds within the galaxy. The highest ionization 
level lines are stronger where collisional ionization dominates. 
Hell 1640, which is a recombination line, increases with both 
the temperature of the gas and the photoionization flux. 

For 0211-122 we group the shock velocities into two pro¬ 
totypes. Vs ~ 1000 km s 'representing high velocity shocks fits 
the UV lines, while Vs < 500 km s 'reproduces the lines in the 
optical range. The results of selected models appear in cols. 4, 
5 and 6. The weighted sum of the line fluxes appears in col. 7 
of Tabled Models MHl and MH2 which refer to Vs = 1000 
km s 'and 500 km s ', respectively, are shock dominated (SD, 
F=0), while model MH3 with Vs =300 kms 'is radiation dom¬ 
inated (RD) H-shock. Model MH2 is included because, overpre¬ 
dicting the CIV and [Oil] 3727 lines, improves the approxi¬ 
mation of the averaged spectrum. O/H is solar and N/H seems 
lower than solar by a factor of ~ 3. We refer to the solar abun¬ 
dances by Allen (1976). They are discussed by Contini (2014a). 
The average model overpredicts the OIV]/Ha , NIV]/Ha and 
CIII]/Ha multiplet ratios to Ha because the calculations account 
for the sum of all the multiplet terms. We do not know which 
terms are included in the observation data. Different geometri¬ 
cal thickness of the clouds D indicate fragmentation by shock 
in a turbulent regime. The 0406-244 spectrum is satisfactorily 
reproduced by the weighted sum (av in column 11) of models 
MH4, dominated by the AGN radiation, and MH5, dominated by 
the shock. The relative weights were calculated phenomenolog¬ 
ically in order to obtain the best approximation to the spectrum, 
even if some line ratios (e.g. [NeV]) are less fitted. Also for this 
galaxy a SD model (MH5) corresponds to the highest Vs (700 


kms^'). More spectra of selected galaxies from the Humphrey 
et al sample are shown and modelled in Table 3. 

The FWHM of the line profiles range between 600 and 900 
kms“'. They are not reported in the table for sake of space. 
In Table 3 we refer all the spectra to Hy8 =1 in order to eas¬ 
ily compare observed line ratios to model results from differ¬ 
ent galaxies. The spectrum presented for 0529-549 in Table 
3 shows a few line ratios. Most are upper limits. Both [OIII] 
5007-I-/H/1 line ratios calculated by high Vs and low Vs mod¬ 
els reproduce satisfactorily the data. [NII]/H/1 and [SII]/HjS de¬ 
pend strongly on the N/H and S/H relative abundances, respec¬ 
tively. So, an averaged model is not significant. The modelling 
of 0828 suggests that carbon should be depleted by a factor of 
~ 5 from the gaseous phase, most probably trapped into grains 
throughout the nebula represented by model MH8. Neon can¬ 
not be included into dust grains because of its atomic structure, 
then we determine the relative weights of the averaged model 
av by fitting the [NeV]/Hy6 line ratio. Model av (col. 8) shows 
that, in contrast to our previous hypothesis, C/H solar is suit¬ 
able to the high velocity gas component of 0838. The observed 
lines of 1138-262 are few and the oxygen lines are well approx¬ 
imated by both the AGN and SB dominated models, therefore 
the relative weights in an eventual average spectrum are less 
constrained. Galaxy 4C-1-23.56 spectrum shows a non negligi¬ 
ble [NeV], while model MHO yields [NeV]=0, therefore the 
SB dominated model MHO is more appropriated. 

The main conclusions by Humphrey et al. are in agreement 
with our modelling. Namely, 1) AGN photoionization with U 
varying between the objects is generally indicated. However, we 
find an SB contribution to 1138-262 and 4C-1-23.56 spectra; 2) 
single-slab photoionization models are unable to reproduce sat¬ 
isfactorily the high and low ionization lines simultaneusly; 3) 
shocks alone do not provide a valid explanation but must be 
added to a photoionization source; 4) N/H is nearly solar for all 
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Table 5. Modelling the spectra from the observations of Vernet et al (2001) 


Flux* 

4C+03.24 

MVl 

0943-242 

MV2 

MV3 

0731+438 

MV4 

MV5 

av 

4C-00.64 

MV6 

z 

SFR2 

3.560 

192.3 


2.922 

100.6 



2.429 

128.6 




2.360 

73.27 


OVI 

- 

- 

5.6 

4.3 

9.5 

- 

- 

- 

- 

- 

- 

ov 

4.8+0.8 

7.3 

- 

- 

- 

- 

+ 

- 

- 

- 

- 

Lya 

254.4±7.3 

237. 

213.8+1.8 

240. 

263. 

428.1+2.2 

426. 

455. 

445. 

261.6+3.4 

233. 

NV 

4.8+0.8 

5.2 

11.2+0.6 

4.4 

14.2 

7.1+0.6 

8.0 

7.7 

7.8 

12.4+1.2 

15.3 

CII 

4.6+1.3 

1.6 

2.4 

0.7 

2.2 

2.2 

1.4 

0.7 

0.94 

1.1 

1.4 

OlV+SilV 

6.1+2.5 

5.0 

8.9 

9.1 

16. 

9.7 

13.0 

3.9 

7.0 

10.4 

14. 

NIV] 

- 

- 

2.4 

4. 

6.8 

1.8+0.3 

3.0 

2.5 

2.67 

- 

- 

Sill 

- 

- 

1.2 

- 

- 

0.9 

- 

- 

- 

- 

- 

CIV 

12.5+2.2 

10.5 

46.5+1.1 

55.2 

40. 

46.5+1.0 

17. 

47.3 

36.9 

29.1+2 

30.07 

[NeV] 

- 

- 

- 

0.4 

0.62 

- 

- 

- 

- 

- 

- 

Hell 

5.2+1.3 

6.0 

47.1 + 1.3 

71.8 

22.8 

30.4+0.5 

23. 

19.6 

20.8 

21.5+0.5 

11.2 

OIII] 

- 

- 

7.6 

4. 

6. 

6.4 

1.7 

2.4 

2.16 

3.6 

4. 

NIII] 

- 

- 

1.4 

2. 

2.9 

- 

0.5 

1.2 

0.96 

- 

- 

Sill 

- 

- 

1.7 

- 

- 

1.4 

- 

- 

- 

- 

- 

Silll] 

- 

- 

5.4 

1.4 

1.6 

2.7 

0.1 

1.1 

0.8 

- 

- 

CIII] 

- 

- 

29.3+3.1 

24. 

14. 

21.2+0.5 

1.4 

24.8 

16.9 

9.6+0.7 

11.7 

Ho- 

30.5 

- 

25.66 

- 

- 

51.36 

- 

- 

- 

31.4 

- 

Vs 3 

- 

130 

- 

100 

130 

- 

140 

140 

- 

- 

120 


- 

60 

- 

100 

150 

- 

150 

130 

- 

- 

60 

F5 

- 

- 

- 

2.5 

- 

- 

- 

6. 

- 

- 

- 

Y 6 

- 

4.55 

- 

- 

6 

- 

5 

- 

- 

- 

4.5 

u 

- 

0.7 

- 

- 

0.15 

- 

0.9 

- 

- 

- 

0.7 

d’’ 

- 

0.013 

- 

0.01 

0.0008 

- 

0.001 

0.1 

- 

- 

0.01 

C/H** 

- 

3.3 

- 

3.3 

3.3 

- 

3.3 

3.3 

- 

- 

3.3 

N/H** 

- 

1. 

- 

1. 

1. 

- 

1. 

1. 

- 

- 

1. 

0/H** 

- 

6.6 

- 

6.6 

6.6 

- 

6.6 

6.6 

- 

- 

6.6 

HyS® 

- 

0.009 

- 

0.012 

0.009 

- 

0.014 

0.13 

- 

- 

2.6e-4 

w 

- 

- 

- 

- 

- 

- 

0.97 

0.03 

- 

- 

- 

ergcm“^ s“^; 

^ in Mo yr“*; 

^ in km s“ 

14* -3 

; in cm ; 

^ in 10^^ photons cm ^ s ^ eV ^ 
calculated in ergcm“^ 

at the Lyman limit; 

in 10^ 

K; ■* in 10*’ 

cm; * in 10 


units; ® 


Table 6. Modelling the spectrum of the SST24 J1434110+331733 nebula (Dey et al 2005), CDF-S202 (Norman et al 2002) and 
Lynx arc (Fosbury et al 2003) 


flux 

SST24* 

MDR 

MDOS 

MDSB 

av 

CDF-S202’ 

FWHM* 

MNR 

MNS 

Lynx'* 

MFO 

MFSB 

z 

2.656 





3.7 




3.36 



SFR(Mo yr-‘) 

134.6 





136.2 




177.8 



Lycr 1216 

28.9+0.15 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Lycr 1216 

4.05+0.01 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Lya 1216 

3.10+0.01 

3.03 

3.3 

4.9 

3.39 

16.4 

<1130 

28. 

22.6 

10.8 

31.8 

38. 

NV 1240 

<0.05 

0.017 

0.03 

0.1 

0.03 

5.9 

<1680 

0. 

8.75 

<0.09 

0.3 

0.001 

SilV 1397 

- 

- 

- 

- 

- 

- 

- 

- 

- 

<0.09 

0.09 

- 

NIV]1483+ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.42 

0.15 

0.05 

CIV 1549+ 

0.417+0.001 

0.6 

0.47 

0.4 

0.49 

9.9 

<1680 

3.19 

12.2 

3.65 

3.6 

0.46 

Hell 1640 

0.4+0.004 

0.82 

0.39 

0.3 

0.45 

2.8 

<680 

0.7 

6.48 

0.11 

0.11 

0.12 

OIII] 1661+ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.56 

0.48 

0.09 

NIII] 1750 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.18 

0.1 

0.12 

Silll] 1883+ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

0.15 

0.14 

0.06 

CIII] 1909 

0.05+0.01 

0.0013 

0.18 

0.03 

0.13 

- 

- 

12.8 

4.8 

0.59 

1.4 

1.0 

[Oil] 3727+ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

<0.25 

0.17 

0.57 

[Nelll] 3869+ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

<0.91 

0.2 

1.16 

Hell 4686 

- 

- 

- 

- 

- 

- 

- 

- 

- 

<0.22 

0.01 

0.02 

W 

- 

0.1 

0.1 

0.1 

0.1 

1 

- 

1 

1 

1 

1. 

1 

[OIII] 5007+ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

10.1 

0.7 

9. 

5084 

0.52+0.004 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Ha 6563 

0.371 

0.3 

0.3 

0.3 

0.3 

1.97 

- 

3. 

3. 

- 

3.7 

3.1 

Vs (kms^*) 

- 

1700 

1700 

300 

- 

- 

- 

1600 

1600 

- 

320 

100 

no (cm“’) 

- 

110 

no 

150 

- 

- 

- 

120 

120 

- 

2000 

300 

F5 

- 

100 

- 

- 

- 

- 

- 

1 

- 

- 

- 

- 

'p 6 

- 

- 

- 

2.5 

- 

- 

- 

- 

- 

- 

- 

1.1 

u 

- 

- 

- 

8 

- 

- 

- 

- 

- 

- 

- 

0.3 

D (10*’ cm) 

- 

300 

3. 

5. 

- 

- 

- 

100 

15 

- 

1 

4900 

c/h’ 

- 

3.3 

3.3 

3.3 

- 

- 

- 

3.3 

3.3 

- 

3.3 

3.3 

n/h’ 

- 

1. 

0.1 

1. 

- 

- 

- 

1. 

1. 

- 

0.2 

1.0 

0/H’ 

- 

6.6 

6.6 

6.6 

- 

- 

- 

6.6 

6.6 

- 

6.6 

3.6 

up calc* 

- 

0.0336 

0.024 

5.68 

- 

- 

- 

200. 

0.0085 

- 

0.15 

8.7 

W 

- 

0.137 

0.822 

0.041 

- 

- 

- 

- 

- 

- 

- 

- 


’ 10 ** ergcm *(Dey et al.); ^ 10 **ergcm '; (Norman et al (2002); ^ in km s *;^inl0 erg cm '(Fosbury et al 2003); ^ in 10*® photons cm ^s ’ 

eV“' at the Lyman limit; ® in 10^ K ^ in 10“^ units; * in erg cm“^ s“' 
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objects and does not vary as much as by a factor of 2-3; 5) The 
line FWHM and the ionization state are the result of interaction 
between the clouds in terms of shocks, when the radiation from 
the AGN or the SB is absent. 


3.2. HDFX 28, a spiral galaxy at z=2.011 from the 
observations of Dawson et al. (2003) 

In search for type 2 QSRs, HDFX 28 was identified as a face-on 
spiral galaxy at z=2.011. Initially, it was considered as an ex¬ 
tended radio source powered by star formation (Richards 2000), 
but the hard X-ray source and the spectroscopy demonstrated 
that this galaxy contains an obscured type 2 AGN (Dawson et al 
and references therein). The galaxy was observed by ISOCAM 
and Chandra. Table |4] shows that a SB dominated model (MDS) 
roughly reproduces the data. Dawson et al spectrum contains 
too few narrow lines, [Nil] 6548, 6548 and Ha to provide a 
consistent modelling of the low Vj gas. On the basis of [Nil] 
6548-i-6584/Ha , we can predict that, if the reddening is negli¬ 
gible and Ha /lij3 -3, [NII]/Hj0 = 4.11, which indicates for a 
SB model U<0.1, T. = 1 -4 1 O^K, =220 km s^', no=200 cm^^, 
D=10'^ cm and solar abundances. For an AGN model, such line 
ratios can be obtained with N/H solar. Vs =100-250 kms^'and 
F< 10^ photons cm“^ s“^ eV ' at the Lyman limit , and for 
Vs =500 kms^^and F < 10'° photons cm"^ s“^ eV“' at the 
Lyman limit. Such [NII]/Hj8 can be also obtained by a SD (F =0) 
model with Vs =300 km s 'The preshock densities are 100-300 
cm“^ and the geometrical thickness of the emitting clouds are 
between 10'^ and 10*^ cm. In conclusion, the [NII]/HjS line ratio 
alone cannot constrain the narrow line model. 

In Fig. 2 we present the modelling of the continuum SED. 
The data are given by Dawson et al 2003, table 1). The contin¬ 
uum SED calculated by the code SUMA shows two curves for 
each model. One represents the bremsstrahlung from the gas and 
the second one, in the IR, represents reprocessed radiation from 
dust. The data in the radio - UV - soft-X-ray range are well fitted 
by the high Vs, SB dominated model MDS. Harder X-ray data in 
Eig. 2 could be reproduced by a model calculated by a shock ve¬ 
locity even higher than 1500 km s“', as suggested by the EWHM 
of some line profiles. The AGN and the SB dominated models 
calculated by Vs =220 km s ', which were predicted on the basis 
of the [NII]/HjS line ratio, reproduce the data in the optical range, 
in agreement with the Ha narrow line EWHM. The modelling of 
the HDE galaxy continuum SED constrains the choice of the low 
velocity model, but it cannot select between an AGN or a SB. 

3.3. Spectrum from CXO 52 by Stern et al. (2002) 
observations 

CXO 52 was identified as an optical, colour selected, high red- 
shift source (Stern et al 2002c) and by optical follow up of X-ray 
sources in the Lynx field (e.g. Stanford et al 1997) . Stern et al 
(2002) report on observations of a type 2 quasar at z=3.288 that 
was identified as a hard X-ray source by Chandra X-ray obser¬ 
vatory. The optical spectrum of CXO 52 was obtained by Keck 
I telescope. Line profile EWHM are broad (> 1000 km s '). A 
shock dominated model with Vs =1500 kms 'and pre-shock 
density no=110 cm“^, resembling those found in other galax¬ 
ies at similar redshifts, fits most of the UV lines (MSS Table 4). 
The Lya /H/? is overpredicted by a factor of 1.6. OVI1035/Hy6 is 
underpredicted and NV 1240/H;8 is well reproduced. The max¬ 
imum temperature of the gas downstream, near the shock-front 
depends on the shock velocity (Sect. 2) therefore, MSS model 
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Fig. 2. Modelling the continuum SED observations Dawson et 
al. (2002, table 1 and references therein). Blue lines:the shock 
dominated model MDS model; black lines : the AGN radiation 
dominated model calculated by Vs =220 km s^', no=200 cm“^, 
logF'=9, and solar abundances; red lines : the SB radiation dom¬ 
inated model calculated by T» =3 lO^K and t/=0.L 
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fits NV/HyS because NV is a high ionization level line. OVI 
which shows a EWHM =2640 km s^'cannot be well reproduced 
by model MSS. 

The Hy6 and the forbidden lines [OIII]5007 h- correspond to 
lower EWHM (170-430 kms^'). The lines are emitted from 
clouds within the galaxy with different Vs . A radiation domi¬ 
nated model (MSR) with Vs =200 kms 'and no=350 cm“^ re¬ 
produces the UV lines (except OVI and NV) within a factor of 2 
and explains the high [OIII]5007 -f4959/Hj 8 . The flux from the 
AGN corresponds to a Seyfert 2. The results of the two models 
are not summed up because they show very different gas condi¬ 
tions in (most probably) well separated regions throughout the 
galaxy. 

To cross check the models we compare the bremsstrahlung 
calculated by models MSS and MSR with the data presented by 
Stern et al (in their table 1) in Eig. 3. The datum at 2-10 keV 
is well reproduced by model MSS, calculated by a high shock 
velocity. The bremsstrahlung maximum depends on the maxi¬ 
mum temperature of the gas downstream. Dust grains can be 
completely destroyed by sputtering at high shock velocities (> 
1000 kms '). The grains in the nebulae represented by models 
MDS and MSS are characterised by a relatively large initial ra¬ 
dius (agr=l /rm) in order to survive to strong sputtering through¬ 
out a high velocity shock-front. However, they can be completely 
destroyed by sputtering and evaporation, as happens for model 
MSS, depending on the cooling rate downstream. The relatively 
low Vs which characterises the gas conditions in model MSR 
yields grain survival. The bremsstrahlung calculated by model 
MSR, which corresponds to gas ionised and heated by the flux 
from the AGN (H-shocks), roughly fits the data in the UV-optical 
and radio ranges. The reprocessed radiation by dust in the IR is 
constrained by only one datum which leads to dig - 10“"' by 
mass. 
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Fig. 3. Modelling the continuum SED observations Stern et al. 
(2002, table 1 and references therein). Blue lines:the shock dom¬ 
inated model MSS model; red lines ; the radiation dominated 
model MSR; 


3.4. The spectrum observed from Q2343-BX418 by Erb et al 
( 2010 ) 

In column 9 of Table 4 we present the observed line ratios for 
BX418 galaxy at z=2.3, defined by Erb et al., as ’’one of the 
youngest and lowest stellar mass continuum-selected z~2 galax¬ 
ies”. 

The spectral lines show different EWHM for different wave¬ 
lengths decreasing from ~ 840 kms 'for Lya to 640 km s 'for 
Hell 1640, down to < 260 km s 'for lines in the optical-near- 
UV range (Table 4, col 10). The observations show that differ¬ 
ent components contribute to different lines. Although permitted 
line fluxes (e.g. Lya and Hell) depend strongly on the photoion- 
izing radiation reaching the nebula, they also increase with the 
temperature of the emitting gas (see Sect. 3.8). Therefore, mod¬ 
els accounting for shock-t photoionization are used. The models 
are presented in Table 4. We first investigate whether an AGN 
(hidden?) is present. The high Lya and Hell 1640 EWHM sug¬ 
gest shock velocities of ~ 850 and ~ 600km s“', respectively. 
They are accompanied by a AGN flux characteristic of Seyfert 2 
galaxies. The densities are high throughout the nebula to reduce 
the primary and secondary fluxes rapidly enough. The temper¬ 
ature drop will lead to very low line intensities in the optical 
range. Actually, broad optical lines were not observed. We could 
not obtain acceptable results with models corresponding to out¬ 
flow of the gas, on the contrary, gas inflow towards the AGN is 
indicated (models MEl-0 and ME2-0). We reproduced the near 
UV and optical line ratios by models corresponding to the NLR 
conditions. Models ME3 and MESB roughly fit the [OII]/Hy6 and 
[OIII]/Hb line ratios. ME3 refers to an AGN, while MESB refers 
to a SB radiation source. Both shows outlow of gas from the 
galaxy and are similar to those found in other galaxies at those 
redshifts. We suggest that all the gas conditions represented by 
the models in Table 4 are present within BX418, contributing 
to the single lines in the spectrum. The models are not summed 
up, in order to show each contribution to each line. BX418 has 
been placed among metal-poor galaxies (Erb et al). However, 
the low metallicity (0.05 solar) was calculated by Erb et al by 
the strong-line metallicity diagnostic as well as by the direct 
methods. They found 12H-log(O/H)=7.9+0.2, while the results 



Fig. 4. The continuum SED modelling of Dey et al. data. Blue 
lines;the AGN model; red lines : the SB model; black lines: the 
shock dominated model and reradiation by dust at 1000 K 


of detailed modelling yield solar 12-i-log(0/H)=8.812. The dis¬ 
crepancy between metallicities calculated by direct methods and 
modelling has been discussed by Contini (2014b). 

3.5. Spectra by Vernet et al (2000) atz~ 2.5 and the 
diagrams 

Vernet et al (2001) present spectropolarimetry of powerful ex- 
tragalactic radio sources at z~2.5 by Keck II telescope, ’’when 
quasars were much more common than now”. Although dust- 
reflected quasar light dominates the UV continuum, they did 
not exclude star formation. Vernet et al. comparing the data 
throughout the NV/Hell vs NV/CIV diagrams by Hamann & 
Eerland (1993) could not easily explain the observed NV/CIV 
and NV/Hell line ratios. In Sect. 4 we will discuss Vernet et 
al and other survey observations throughout the Hell/CIV vs 
Lya /CIV diagrams, where the data are compared with model re¬ 
sults calculated by SUMA. UV line ratios are generally referred 
to CIV which is one of the strongest lines. Moreover, we have 
found that the N/H relative abundances throughout galaxies are 
scattered by a factor of ~ 10 (Contini 2014a, fig. 5), therefore we 
suggest that NV lines are less adapted to diagnostic diagrams. 

We present in Table |5] the results of shockH-photoionization 
models for some of the Vernet et al spectra. 4C-1-40.36 and 
4C-1-48.48 were excluded because the spectra are most probably 
blended by a foreground HII galaxy at z=0.404 and at z=0.684, 
respectively. The modelling results indicate that composite mod¬ 
els referring to the SB-i-shocks, to the AGN-i-shocks should be 
adopted and that different physical conditions from different re¬ 
gions contribute to the spectra. All the models used to reproduce 
Vernet et al data in Table 5 show relatively low Vs and no, which 
are reflected in low SER. Vernet et al found a SER range be¬ 
tween ~ 60 Mq yr“' to < 2 Mq yr“'. They claim that in 0731 -1-438 
the only object where reprocessed radiation from AGN does 
not dominate the continuum (for the other objects starlight con- 
ributes less than half of the UV continuum, see Vernet et al 2001, 
table 5) and assuming E(B-V)=0.1 and a stellar population age 
of 10^ yr, Vernet et al calculate SER between 30 and 120 Mq 
yr“'. Our calculations give SER ~ 128 M© yr^'for 0731-1-438, a 
minimum SER ~ 24.4 for 4C-1-23.56a, etc. (see Table 5) 
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Fig. 5. The continuum SED modelling of SFD-S202. The data 
are adapted from Norman et al. (2002, fig. 6) Blue lines: the 
shock dominated model MNS; red lines : the AGN model MNR. 
black solid line : the X-ray data 


3.6. Line and continuum observations of an extended nebuia 
by Dey et at. (2005) 

Massive ellipticals at high z (>3) are surrounded by large (100- 
200 kpc) luminous (>10"^^ erg s^' ) Lya nebulae and small em¬ 
bedded star forming objects. Investigating (>100 kpc) halos and 
environments of luminous radio galaxies affected by AGNs and 
ejecta of similarly large Lya nebulae, Dey et al (2005) reported 
on observations of a very large, spatially extended (160 kpc) 
Lya emitting nebula at z=2.656. Spectroscopic observations of 
SST24 J1434110-1-331733 were obtained by the Low Resolution 
Imaging Spectrometer (LRIS) on the Keck I telescope. 

The spectroscopic observations reveal very complex re¬ 
gions within the extended nebula, including several continuum 
sources. Dey et al suggest that many active galaxies are forming 
in that region which lead to three main contributions to the con¬ 
tinuum spectrum. Galaxy A shows interstellar absorption lines 
and no Lya emission; the spectrum of the second continuum 
source is quite red and shows weak CIV and CIII] emission at 
z=2.656. Star formation is suggested because the source is close 
to MIPS 24 yum location. The third continuum source lies in the 
centre of the Lya nebula associated to CIV and Hell. The veloc¬ 
ity structure of the nebula deduced from the line profiles shows 
that the velocity varies uniformly across the central region. If 
due to gas kinematics, the variation can be due to infall, outflow, 
or rotation. 

We refer to the line and continuum observations by Dey et 
al. The results of modelling are presented in Table 6. On the 
basis of the observed lines corresponding to broad profiles, we 
have run models with Vs > 1000 km s“^ We find a large contri¬ 
butions to the line fluxes from high velocity (Vs =1700 kms^') 
shock dominated gas (models MDOS) and photoionized by an 
AGN (model MDR). Moreover, the contribution from gas with 
Vs =300 kms 'and photoionised by a SB is evident (MDSB). 
The high velocities are at the limit of those observed in AGN 
NLR. The AGN flux is characteristic of Seyfert 2. The starburst 
temperature which leads to a satisfactory fit of the line ratios 
(2.5 10^ K) is as high as that found in galaxies showing activ¬ 
ity (Contini 2014b). Also the ionization parameter is high. The 
gas throughout the SB shows ’’normal physical conditions” in 


Fig. 6. The profiles of the physical quantities throughout the 
emitting nebula for the SD model MLO. The shock-front is 
on the left. Top : electron temperature (solid black), electron 
density (dashed black). Middle : H'^/H (solid black), He'^^/He 
(dashed black), C^^/C (dashed red), C^^/C (solid red), N^+ZN 
(solid blue), 0^/0 (dashed green), O^^/O (solid green). Bottom 
: log (Lya /Hy6 ) (solid black), log(Ha jWfi ) (dashed black) and 
log(Lya /Ha ) (dot-dashed red). 

the emitting nebulae of galaxies at such redshifts. However, to 
have a good fit of the NV 1240 line, we have adopted N/H lower 
than solar by a factor of ~ 10 in the shock dominated gas. This 
indicates that it is located in a different region of the galaxy. The 
weighted sum of the calculated spectra (av) appears in Table 6 
col. 6. We suggest that the A 5084 line refers to scandium ("F 3 / 2 - 
> White 1934, fig. 14.7). Scandium can be present is SN 
ejecta (e.g. Bose et al. 2015). 

Fig. 4 shows that the continuum SED calculated by the 
models selected by fitting the line ratios reproduce the ob¬ 
served continuum data presented by Dey et al. (2005, fig.2). 
Bremsstrahlung from gas heated by the strong shock shows a 
maximum at high frequency (Contini, Viegas, Prieto 2004). X- 
ray emission is predicted. The mass of the emitting gas (6 10*^ 
Mq ) was calculated by Dey et al. The mass of dust can be cal¬ 
culated by dig - 0.003 by mass which results from the fit of the 
continuum SED in the IR. 

3.7. A ctassic type 2 QSO at z=3.7 ; CDF-S202 observations 
by Norman et at (2002) 

Norman et al observed the AGN detected at z=3.7 by the 
Chandra Deep Field-South 1 Ms exposure. It is a distant type 
2 AGN with UV showing FWHM of ~ 1500 km s *. The optical 
spectrum was obtained with the multislit mode on the FORSl 
on the ESO/VLT-ANTU. Detailed modelling (Table 6, cols. 9, 
10) indicates that a shock dominated model (MNS) best fits 
the lines within a maximum discrepancy factor of 2.3 for Hell 
1640 (Table 6 last column). This model calculated by Vs =1600 
km s 'would provide OVI 1035/14/1 ~ 200. Norman et al found a 
high metallicity, while we suggest N/H by a factor of 1.5 lower 
than solar and C/H about solar, considering the uncertainty of 
the line fluxes. Referring the line to CIV=10 instead of H/1, as 
it is generally done for UV spectra, we would obtain an accept¬ 
able fit of all the lines. The radiation (AGN-i-shock) dominated 
model (MNR) calculated by the high Vs reproduces the line ra- 
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log distance from shock front [cm] 


Fig. 7. Same as in Fig. 6 for model MFSB 



log distance from the opposite edge[cm] 


tios roughly. In Fig. 5 we present the modelling of the contin¬ 
uum SED. The data were adapted from Norman et al. fig. 6. 
The continuum calculated by the radiation (-Hshock) dominated 
model MNR reproduces the data in the optical- UV better than 
the shock dominated model MNS, so the MNR model should be 
accounted for, but with a low relative weight. 

3.8. Spectrum of the Lynx arc at z=3.357 by Fosbury et al 
(2003) 

The Lynx arc (z=3.357)is characterised by a very red R-K colour 
and strong, narrow emission lines. Fosbury et al claim that anal¬ 
ysis by HST WFPC2 imaging and Keck optical and IR spec¬ 
troscopy indicate a HII galaxy magnified by a factor of 10 by a 
cluster environment including stars with T, ~ 10®K, U ~ I and 
low metallicity (0.05 solar). In Table 6, col 11 the observed spec¬ 
trum of the Lynx arc presented by Losbury et al (2003, table 2) is 
reported. The spectrum shows a relatively large number of line 
ratios including [OII]/HyS and [OIII]/Hy6 in the optical range, re¬ 
produced by Losbury et al by a pure photoionization model with 
T* =8 10"^ K logt/=-l and Z=0.05Zo, which was suggested by 
comparison with synthetic SEDs calculated by STARBURST99 
(Leitherer et al 1999). The LWHM of the profiles are < 100 
km s 'for lines in the optical and higher for those in the UV. We 
searched by modelling the spectrum for a (hidden) AGN, but 
AGN dominated models lead to HeII/Hj6 > 1, higher by a factor 
>10 than observed. The shock dominated model MLO (Table 
6, col. 12) calculated by a relatively high shock velocity (320 
km s“') and high preshock density (no=2000 cm“^) fits most of 
the UV line ratios but underpredicts [OIII] 5007-i-/HjS by a fac¬ 
tor > 10. A SB radiation dominated model MLSB, calculated by 
T. =1.1 lO^K, t/=0.3,D=4.9 10^° cm, V^ =100 km s^', no=300 
cm“^ and O/H ~ 0.5 solar, reproduces satisfactorily most of the 
line ratios from relatively low ionization levels. Model MLSB 
which accounts for photoionization by a SB -H shocks represents 
gas ejected outwards, while for model MLO the direction of the 
emitting gas is ambiguous. 

Linally, we show in Ligs. 6 and 7 the profiles of the physical 
parameters throughout the emitting nebula for models MLO and 
MLSB, respectively, in order to understand the results. In Ligs. 6 
and 7 the electron temperature and electron density appear in the 
top diagrams. In Lig. 7 the nebula is devided into two halves in 
order to obtain a comparable view at the two edges. Logarithmic 


symmetric scales are used for the X-axis. The shock front is on 
the left of the left panel and the right edge of the right panel 
is reached by the radiation flux from the SB. In the middle di¬ 
agrams of Ligs. 6 and 7 the fractional abundances of the ions 
corresponding to significant lines are shown. The bottom dia¬ 
grams are dedicated to the Lya fli/3 , Ha /H/S and Lya /Ha line 
ratios (in logarithm). Ligs 6 and 7 show that Lya /Ha follows 
the H^ recombination trend. The two sides of the nebula (shock 
dominated at the left and radiation dominated at the right of Lig. 
7) are bridged by diffuse secondary radiation. In fact, the tem¬ 
perature is ~ IO^'K even in the extended internal region of the 
nebula. Most of the C, N, O lines come from the radiation dom¬ 
inated side of the nebula. 


3.9. Finkelstein et al survey 

Linkelstein et al (2011) present in their table 1 the measured line 
fluxes for a sample of galaxies at 0.24 < z < 0.324 by spec¬ 
troscopy of LAE using the Hectospec spectrograph with the 6 m 
MMT telescope. 

We report in Table 7 the results of modelling the observed 
line ratios and we add the UV lines NV, CIV and He II consis¬ 
tently calculated by the same models. These lines which are not 
observed by Linkelstein et al will be discussed in the following. 
In Table 7 the observed line ratios (EGS7-EGS25, corrected by 
Ha /HyS =3) are followed in the next rows by the results of SB 
models (MSB7-MSB25)) and AGN models (MAG7-MAG25)). 
The input parameters are reported in Tables 8 and 9). 

Lya lines are calculated consistently with the other lines. 
Lya /Hy8 calculated line ratios (col. 3) reproduce the data within 
a factor of 2. In col. 4 the SLR are indicated for each galaxy. 
The percentage uncertainty appears in col. 5. The galaxies in the 
Linkelstein et al sample present rather homogeneous character¬ 
istics. The O/H relative abundances are at most solar while the 
N/H ratios are lower than solar by a factor <10. Densities and 
velocities are relatively low. The flux from the active nucleus is 
in the range of that found in low luminosity AGNs. The tem¬ 
perature of the stars does not show the maxima corresponding 
to outbursts (Contini 2014b). The geometrical thickness of the 
clouds shows that fragmentation is present in the SB neighbour¬ 
hood rather than in AGNs. 
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Table 7. Modelling the corrected line ratios to =1 observed by Finkelstein et al (2011) 


ID 

z 

SFR 
Mo yr-‘ 

Lya 

1215 

% 

NV 

1240 

CIV 

1550 

Hell 

1640 

[Oil] 

3727+ 

[OIII] 

5007+ 

Ha 

6563 

[Nil] 

6584 

EGS7 

0.2440 

0.4 

17.53 

15.9 

- 

- 

- 

2.14 

5.65 

3. 

0.134 

MSB7 



24.3 

- 

0.0014 

0.34 

0.078 

2.3 

5.6 

2.93 

0.16 

MAG7 



28.6 

- 

0.0009 

0.15 

1.22 

2. 

6.4 

3.3 

0.18 

EGS8* 

0.2395 

0.035 

16.7 

32.8 

- 

- 

- 

2.47 

3.7 

3. 

<0.03 

MSB8 



23.8 

- 

0.0024 

0.28 

0.0037 

2. 

3.5 

2.95 

0.1 

MAG8 



28.6 

- 

0.0005 

0.14 

0.89 

2. 

3.46 

3.34 

0.1 

EGSIO 

0.2466 

0.12 

33.9 

16.7 

- 

- 

- 

3.06 

4.37 

3. 

<0.06 

MSB 10 



24.6 

- 

0.016 

0.46 

0.044 

3.2 

4.5 

2.94 

0.12 

MAGIO 



28. 

- 

0.001 

0.029 

1.1 

3.1 

4.65 

3.2 

0.14 

EGSll 

0.2524 

0.059 

41.1 

35.3 

- 

- 

- 

2.78 

1.53 

3. 

0.84 

MSB 11 



24.5 

- 

0.1 

0.3 

0.0026 

2.3 

1.5 

2.96 

0.87 

MAGll 



26.5 

- 

0.0027 

0.1 

0.67 

2.5 

1.57 

3.3 

0.9 

EGS12 

0.2515 

0.062 

48.7 

20.8 

- 

- 

- 

2.8 

2.7 

3. 

0.38 

MSB 12 



24.5 

- 

0.067 

0.42 

0.01 

2.9 

2.63 

3. 

0.46 

MAG12 



28. 

- 

0.0027 

0.28 

0.8 

2.5 

2.8 

3.3 

0.37 

EGS13 

0.2607 

0.46 

17.96 

21. 

- 

- 

- 

2.8 

2.26 

3. 

0.43 

MSB 13 



24.5 

- 

0.067 

0.42 

0.0092 

2.9 

2.45 

2.95 

0.47 

MAG13 



31.54 

- 

0.017 

0.25 

0.73 

3. 

2.35 

3.3 

0.46 

EGS19 

0.2666 

0.11 

33.96 

21.9 

- 

- 

- 

2.4 

1.23 

3. 

0.86 

MSB 19 



24.4 

- 

0.1 

0.31 

0.0046 

2.36 

1.3 

3. 

0.86 

MAG19 



26.4 

- 

0.0027 

0.1 

0.61 

2.5 

1.25 

3. 

0.96 

EGS20* 

0.2680 

0.073 

11.57 

27.2 

- 

- 

- 

3.04 

4.2 

2.9 

<0.137 

MSB20 



24.3 

- 

0.024 

0.46 

0.015 

2.96 

4.12 

2.93 

0.14 

MAG20 



28. 

- 

0.001 

0.29 

1.1 

3. 

4.5 

3. 

0.14 

EGS21 

0.2814 

0.06 

46.4 

21.6 

- 

- 

- 

1.98 

0.415 

3. 

0.12-1.4 

MSB21 



32. 

- 

0.063 

0.155 

0.0015 

2.2-1.8 

0.47-0.4 

3.3 

0.13-1.2 

MAG21a 



30.5 

- 

0.002 

0.3 

0.27 

1.8 

0.45 

3.46 

1.5 

MAG21b 



28.3 

- 

0.0027 

0.14 

0.31 

2.1 

0.42 

3.4 

0.15 

EGS22 

0.2828 

0.12 

242. 

35.9 

- 

- 

- 

3.17 

0.74 

3. 

<0.14 

MSB22 



25.5 

- 

0.023 

0.9 

0.011 

3. 

0.75 

2.96 

0.13 

MAG22 



27.9 

- 

0.002 

0.22 

0.5 

3.2 

0.72 

3.3 

0.12 

EGS23 

0.2865 

0.42 

18.4 

15.8 

- 

- 

- 

2.65 

3.7 

3. 

0.32 

MSB23 



24.2 

- 

0.045 

0.55 

0.012 

2.65 

3.74 

2.93 

0.31 

MAG23 



27.1 

- 

0.004 

0.3 

1.0 

2.60 

3.88 

3.2 

0.37 

EGS25 

0.3243 

1.06 

6.5 

16.6 

- 

- 

- 

1.76 

5.94 

3. 

0.16 

MSB25 



24. 

- 

0.04 

0.68 

0.023 

1.8 

6.2 

2.93 

0.15 

MAG25 



28.9 

- 

0.0005 

0.11 

1.13 

1.8 

5.95 

3.3 

0.14 


* not Balmer corrected because the observed Ha /Hj8 ~ 3. 


Table 8. Input parameters for Finkelstein et al. MSB models Table 9. Input parameters for Finkelstein et al. MAG models 


model 

Vs 

no 

T. 

U 

N/H 

O/H 

D 

W 

model 

z 

Vs 

no 

F 

N/H 

O/H 

D 

w 


1 

2 

3 

4 

5 

6 

1 

8 



1 

2 

3 

4 

5 

6 

1 

MSB7 

100 

100 

6.4 

0.02 

0.1 

5. 

9. 

0.02 

MAG7 

0.2440 

100 

150 

7. 

0.2 

5.6 

3 

0.066 

MSB8 

100 

100 

4.2 

0.03 

0.2 

6.6 

3. 

0.0196 

MAG8 

0.2395 

100 

150 

3. 

0.1 

5.6 

3 

0.064 

MSB 10 

150 

100 

5.9 

0.02 

0.1 

6.0 

4. 

0.019 

MAGIO 

0.2466 

100 

150 

3. 

0.1 

6.6 

3 

0.029 

MSB 11 

150 

100 

4. 

0.03 

0.9 

6.6 

4. 

0.027 

MAGll 

0.2524 

100 

150 

1.4 

0.8 

6.6 

3 

0.072 

MSB 12 

180 

100 

4.8 

0.025 

0.4 

6.6 

5. 

0.024 

MAG 12 

0.2515 

100 

150 

2.5 

0.3 

6.0 

1.4 

0.03 

MSB 13 

180 

100 

4.7 

0.025 

0.4 

6.6 

5. 

0.025 

MAG 13 

0.2607 

120 

150 

2.6 

0.3 

6.0 

1.9 

0.045 

MSB 19 

210 

100 

3.9 

0.04 

0.8 

6.6 

5. 

0.038 

MAG19 

0.2666 

100 

150 

1.2 

0.8 

6.6 

3 

0.072 

MSB20 

220 

100 

4.8 

0.035 

0.13 

6.6 

1. 

0.029 

MAG20 

0.2680 

100 

150 

2.9 

0.1 

6.6 

1 

0.03 

MSB21 

170 

100 

4. 

0.015 

0.1 

6.6 

300 

0.058 

MAG21a 

0.2814 

100 

150 

0.77 

0.2 

6.6 

1.6 

0.026 

MSB22 

200 

150 

4.2 

0.015 

0.06 

5.2 

80 

0.023 

MAG21b 

0.2814 

110 

150 

0.8 

0.15 

6.6 

4 

0.068 

MSB23 

190 

100 

4.4 

0.034 

0.2 

5.6 

1. 

0.022 

MAG22 

0.2828 

110 

150 

0.86 

0.08 

6.0 

4 

0.045 

MSB25 

170 

100 

4.9 

0.034 

0.15 

5. 

1. 

0.023 

MAG23 

0.2865 

100 

100 

1.8 

0.3 

6.6 

2 

0.018 

l:in km 

; 2:in cm' 

; 3: 

:in lO^'K 

; 4: - ; 

5:in 10^^ 

; 6:in 

lO-'* ; 

7:in 10*^ cm ; 

MAG25 

0.3243 

100 

150 

6. 

0.15 

5.8 

4 

0.086 


8:in erg cm ^ s ' l:in kms *; 2:in cm 3: in 10^ photons cm ^ s * eV * at the Lyman limit ; 

4:in 10“^^; 5:in 10“^; 6:in lO'^cm ; 7:in ergcm^^ s“' 


3.10. Predicted Lya and UV lines from galaxies at 
0.06<z<0.9 

The R[o///] ([OIII]5007+/[OIII]4363) line ratios indicate that gas 
densities and temperatures in ranges larger than those deduced 
in average from the observations should be accounted for. This is 
valid for the AGN NLR spectra (see e.g. (Contini & Aldrovandi 
1986, fig. 4 and references therein, etc.) and for most of the 
objects throughout the galactic medium (such as SNR, novae. 


symbiotic stars) and the extra-galactic one. In Seyfert type 1 
galaxies [GUI] lines from the BLR are collisionally deexcited. In 
Seyfert type 2 the [OIII] lines correspond to gas at Vs =300-500 
km s^'and densities of 10“*-10^ cm“^. In intermediate type galax¬ 
ies R[o///] ranges between high (> 100) and low values (< 20). 
LINERS (Heckman 1980, Ferland & Netzer 1983, Ho et al 1993) 
cannot be considered as a low F case (Contini & Aldrovandi 
1983, fig. 3). Problems can be solved by adding shock wave hy- 
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Table 10. The data observed by Ly et al (2014) 


ID 

zspec 

FWHM 

[Oil] 

W 

[OIII] 

[OIII] 

Ho- 

[Nil] 

log([OII]/H/l) 

log([OIII]5007+/H/3 



kms“^ 

3727+ 

4861 

4363 

5007+ 

6563 

6583 

- 

- 

MMTOl 

MMT02 

0.6380 

0.4327 

220.0 

247.0 

10.58 

18.20 

9.21 

5.76 

2.25 

0.60 

73.32 

25.83 

0.00 

0.00 

0.00 

0.00 

0.184-- 

0.621«;j«/ 

0.886--^ 

0.644®^ 

MMT03 

0.4809 

256.4 

20.40 

5.81 

0.95 

12.23 

0.00 

0.00 

u. /:?j0 2oo 

0 774O.007 
' '^0.008 

MMT04 

0.3933 

248.6 

49.30 

21.70 

1.38 

130.05 

0.00 

0.00 

0 355^-^^^ 

^•■^‘^^0,009 

MMT05 

0.3846 

232.1 

3.70 

7.02 

1.29 

57.92 

0.00 

0.00 

-0 266^- 

^•^^^0.034 

MMT06 

0.3995 

260.2 

21.33 

6.04 

0.83 

21.97 

0.00 

0.00 

0.734!!®® 

0.543j!'“^ 

MMT07 

0.3896 

238.0 

19.84 

20.84 

2.25 

153.87 

0.00 

0.00 

0.102^ 

0.8720™“ 

MMT08 

0.6335 

239.4 

27.33 

10.90 

1.88 

42.30 

0.00 

0.00 

0.613§'i| 

0.584^^^^ 

MMT09 

0.4788 

235.2 

16.74 

6.82 

1.49 

39.48 

0.00 

0.00 

0 5858:?^? 
0.1440«^ 

0 762^:8^^ 
0.824g'gi« 

MMTIO 

0.0683 

286.1 

22.46 

18.19 

2.02 

121.87 

59.04 

1.26 

MMTll 

0.1752 

276.9 

42.42 

22.01 

2.05 

144.89 

78.37 

3.13 

0.386|g 

0.840®™' 

0.8238™| 

0.8128™ 

MMTl 2 

0.6405 

214.5 

23.57 

6.18 

1.31 

40.13 

0.00 

0.00 

MMTl 3 

0.4696 

209.1 

6.03 

6.51 

1.15 

39.59 

0.00 

0.00 

0 072^® 

o.78o8;|;| 

MMT14 

0.4644 

280.6 

100.27 

56.75 

2.01 

354.12 

0.00 

0.00 

u.^OOo,019 

0.7940;<>“ 

Keckl 

0.8390 

124.6 

4.62 

7.80 

0.80 

61.10 

0.00 

0.00 

-0 016° °®® 
La 

0.843;® 

0.676:g 

o-^ss®:®!! 

Keck2 

0.6230 

89.5 

3.49 

2.33 

0.67 

17.60 

0.00 

0.00 

Keck3 

0.7906 

121.2 

27.88 

12.50 

0.27 

59.17 

0.00 

0.00 

Keck4 

0.8829 

178.5 

13.77 

7.28 

0.33 

35.60 

0.00 

0.00 

Keck5 

0.8353 

101.1 

1.85 

1.10 

0.58 

7.99 

0.00 

0.00 

0 362®'^'® 

-loS-^^ 

_'-'”°0.937_ 

0.8790;|5 

__ 

Keck6 

0.8237 

90.5 

0.05 

0.88 

0.19 

5.36 

0.00 

0.00 


Note : the line fluxes are in 10 erg cm ^ s 'with 68 % confidence uncertainties 


log(R[o/;/]) 


1 462'’'"' 

i:549|i 

i.yjUo 141 

i.uj5'o.192 

1 

^•■^^^0.408 
1 7770.132 

i26o8:ii 

,3:^1 

i . J rsQS 


:8:i§8 

'8:^i 

.307 

.112 


1.717 
1.803 
1.318 
1.524 
2.2151 

2i5:1M 

toigSSi 

2.U18d290 

1 

1-3468:11 

'-’^°0.32n 


drodynamics to the calculations of the spectra (Contini 1997 and 
references therein) or by assuming a density gradient or a strati¬ 
fication of emitting filaments (Pequignot 1984, Filippenko 1984, 
etc.) The importance of the [OIII]4363 line recently emerged for 
high z galaxies, because constraining the models. To cross-check 
this issue we calculate the Lya /Hj6 line ratios for galaxies which 
show spectra constrained by the [OIII] 4363 line (see Contini 
2014a, Ly et al.2014). We choose the Ly et al (2014) survey be¬ 
cause the reddening corrected spectra contain [OIII]50073-4959, 
[OIII] 4363, [Oil] 3727, H/3 and Ha lines. 

In Table 10 the results of Ly et al. (2014) observations ob¬ 
tained by the Subaru telescope are reported. We compare the 
calculated with observed line ratios in Table 11. The uncertainty 
for model calculations is about 10 %. 

In Table 12 models ml-m20 are calculated by SB domi¬ 
nated (-Hshock) models, while models mpll-mpl20 (described 
in Table 13) are calculated by AON (-Hshocks) dominated mod¬ 
els. Table 11 show that [OII]/HjS and [OIII]/I^ line ratios are 
well fitted by both SB and AGN models, but, regarding the 
[OIII]5007 h-/[OIII] 4363 line ratios, the SB dominated models 
fit the data of nearly all galaxies, while MMT5, MMT12, and 
Keck3 are better reproduced by AGN dominated models. For 
MMT4, Keck4 and Keck6 the SB and AGN photoionizing mod¬ 
els act with similar weights. 

The [OIII]4363/Hj6 is the key line ratio for choosing between 
SB or AGN dominated models. Indeed the shock plays an impor¬ 
tant role, but the effect of photoionization cannot be neglected. 
The ionization parameter (t/) and the SB effective temperature 
(T* ) for SB and the power-law flux from the AGN are deter¬ 
mined by fitting the [OIII]5007 -h/H/ 3 and [OII]3727/H;S line ra¬ 
tios. With those input parameters, the [OIII]4363/Hj6 line ratios 
for SB are higher by factors very roughly ~ 10 than those calcu¬ 
lated for the AGNs, indicating that the shock effect is relatively 
strong. So [OIII]5007-i-/[OIII]4363 are lower for starbursts. 

The densities in the downstream emitting regions are higher 
by a factor of ~ 10, depending on Vs . Regarding the 
[NII]/I^ line ratios, the Ha [Wfi line ratios for MMTIO and 
MMTll are roughly between 3 and 4, the [NII]/Hj6 line ratios 
were directly calculated. The result of N/H is given only for the 


SB dominated models which well reproduce all the other line 
ratios. We obtained low N/H relative abundances (0.25 10 “^) for 
the two galaxies MMTIO and MMTl 1. The results for O/H (in 
10 4 units) are shown in Tables 12 and 13. O/H relative abun¬ 
dances calculated by detailed modelling are closer to solar and 
higher than those calculated by direct methods. 

4. Discussion and conciuding remarks 

In previous sections we have been investigating the gas physical 
conditions and element abundances in galaxies at relatively high 
z, selected by their characteristic line spectra, i.e. including the 
Lya line and the lines which constrain the models. A detailed 
modelling process has been used. All the models are calculated 
at the emitting nebulae in the different galaxies. They lead to 
Lya /HjS higher than observed, except for some SB dominated 
models, MDS (Table 4) , MVl and MV6 (Table 5) and some 
models in the Finkelstein et al survey (Tables 7-9). The mod¬ 
els best reproducing the data in Tables 2-6 show that most of 
the spectra are fitted by the summed contributions of two types : 
those emitted from the high velocity (~ 700-1500 km s“') clouds 
(generally close to an AGN), which explain the high ionization 
level lines and UV - X-ray data in the continuum SED, and those 
emitted from the relatively low Vj (~ 100-300 kms^'), no (100- 
300 cm“^) radiation dominated clouds (by a AGN or a SB) which 
confirm the narrow optical line profiles, as well as the optical 
spectral and continuum data. In some cases one or more model 
components show that the photoionizing flux from the active 
source is prevented from reaching the gas, or that the source is 
absent. Therefore, the clouds within the galaxies are ionized and 
heated only by the shock. In the HDF galaxy (Dawson et al) the 
high Vs (1300 km s ') clouds are found in the surroundings of a 
SB corresponding to a high T,. In the BX4I8 galaxy (Erb et al), 
the high velocity gas (Vs =600-800 kms^*) shows a preshock 
density higher by a factor >10 than in the other objects. The 
shock velocities were most probably reduced by collision with 
high density gas. These clouds are located in the surroundings 
of an AGN. This is valid also for the shock dominated clouds 
(model MEO, Table 6) in the Lynx arc (Eosbury et al). A few 
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Table 11. Modelling the corrected line ratios to =1 observed by Ly et al. (2014) 


ID 

OVI 

1034 

ov 

1215 

'Lya 

1215 

NV 

1240 

CIV 

1549 

Hell 

1640 

[Oil] 

3727+ 

[OIII] 

4363 

[OIII] 

5007+ 

MMTOl 

- 

- 

- 

- 

- 

- 

1.53 

0.26 

7.69 

ml 

38.7 

7.86 

25.89 

3. 

7.82 

0.26 

1.53 

0.26 

7.69 

mpll 

5.3 

1.11 

32. 

0.43 

0.68 

1.45 

1.5 

0.05 

7.8 

MMT02 

- 

- 

- 

- 

- 

- 

4.18 

0.124 

4.4 

m2 

57.8 

10.6 

26.2 

3.75 

9.7 

0.14 

3.9 

0.23 

4.2 

mpl2 

5. 

1. 

30.9 

0.39 

0.98 

1.44 

4. 

0.035 

4.74 

MMT03 

- 

- 

- 

- 

- 

- 

5.43 

0.18 

2.18 

m3 

54.9 

10.9 

36.7 

4.2 

10.9 

0.117 

5.2 

0.2 

2.5 

mpl3 

10.3 

1.93 

30.8 

0.65 

1.64 

1.15 

5. 

0.04 

2.4 

MMT04 

- 

- 

- 

- 

- 

- 

2.26 

0.07 

5.94 

m4 

46.8 

9. 

26. 

3.65 

9.28 

0.18 

2.7 

0.2 

5.43 

mpl4 

2.7 

0.6 

30.7 

0.32 

0.79 

1.78 

3. 

0.034 

5.8 

MMT05 

- 

- 

- 

- 

- 

- 

0.54 

0.18 

8.33 

m5 

3.16 

0.67 

24.6 

0.37 

0.96 

0.48 

0.9 

0.04 

8.33 

mpl5 

2.42 

0.52 

32.2 

0.36 

0.91 

2.11 

1.2 

0.07 

8.4 

MMT06 

- 

- 

- 

- 

- 

- 

5.42 

0.167 

3.49 

m6 

69.4 

16. 

49.7 

13.4 

30 

0.425 

6. 

0.25 

3.1 

mpl6 

8.3 

1.52 

30.5 

0.51 

1.31 

1.29 

5.1 

0.038 

3.5 

MMT07 

- 

- 

- 

- 

- 

- 

1.26 

0.12 

7.44 

ml 

36.7 

7.74 

26. 

3.45 

8.73 

0.3 

1.35 

0.19 

7.33 

mpl7 

5.3 

1.11 

32. 

0.43 

0.68 

1.65 

1.4 

0.05 

7.8 

MMT08 

- 

- 

- 

- 

- 

- 

4.1 

0.21 

3.84 

m8 

57.5 

10.5 

26. 

3.73 

9.61 

0.14 

4. 

0.226 

3.94 

mpl8 

6.65 

1.3 

30.9 

0.5 

1.25 

1.28 

4.4 

0.036 

3.88 

MMT09 

- 

- 

- 

- 

- 

- 

3.84 

0.265 

5.78 

m9 

59.2 

11.1 

26.2 

4.2 

10.8 

0.18 

3.5 

0.264 

5.9 

mpl9 

4.1 

0.83 

31.2 

0.37 

0.92 

1.66 

3.71 

0.04 

5.75 

MMTIO 

- 

- 

- 

- 

- 

- 

1.39 

0.128 

6.67 

mlO 

26.8 

6.35 

25.9 

0.86 

8.75 

0.23 

1.41 

0.16 

6.42 

mpllO 

7. 

1.38 

32.3 

0.54 

0.58 

1.54 

1.3 

0.04 

6.5 

MMTll 

- 

- 

- 

- 

- 

- 

2.43 

0.1 

6.65 

mil 

30.5 

6.63 

25.6 

0.63 

6.75 

0.123 

2. 

0.167 

6.0 

mpll 1 

2.66 

0.52 

32.8 

0.21 

0.53 

1.43 

2.3 

0.04 

6.3 

MMT12 

- 

- 

- 

- 

- 

- 

6.9 

0.31 

6.5 

ml2 

118.7 

21.8 

27.8 

7.7 

19.8 

0.27 

6.7 

0.45 

6.2 

mpll 2 

57.6 

20.4 

27.1 

3. 

8.2 

1.17 

6.2 

0.28 

6.4 

MMT13 

- 

- 

- 

- 

- 

- 

1.18 

0.18 

6.03 

ml3 

26. 

6,25 

25.9 

0.9 

9.12 

0.21 

1.39 

0.15 

6.0 

mpll 3 

4.74 

1.04 

32.3 

0.47 

1.23 

1.53 

1.2 

0.055 

6.28 

MMT14 

- 

- 

- 

- 

- 

- 

1.85 

0.04 

6.22 

ml4 

72.4 

17. 

28.3 

3.12 

30.3 

0.38 

2.4 

0.2 

6.1 

mpll 4 

4.64 

1.17 

26.45 

1.47 

3.3 

2.55 

1.8 

0.059 

6.1 

Keckl 

- 

- 

- 

- 

- 

- 

0.96 

0.12 

7.89 

ml5 

28.9 

6.3 

25.8 

2.8 

7.25 

0.47 

0.9 

0.16 

7.8 

mpll 5 

2.4 

0.52 

32.4 

0.36 

0.9 

2.02 

1.2 

0.063 

7.88 

Keck2 

- 

- 

- 

- 

- 

- 

1.39 

0.25 

6.97 

ml6 

0. 

0.02 

24.3 

0.043 

10.9 

0.165 

1.3 

0.29 

6.7 

mpll 6 

1.41 

0.33 

34.9 

0.23 

0.59 

1.96 

1.27 

0.05 

7. 

Keck3 

- 

- 

- 

- 

- 

- 

4.55 

0.034 

4.74 

ml7 

13.6 

10.4 

27.4 

4.12 

12.9 

0.07 

4. 

0.32 

4.57 

mpll 7 

0.081 

0.3 

26.4 

0.42 

1.81 

1.55 

4.4 

0.049 

4.75 

Keck4 

- 

- 

- 

- 

- 

- 

1.88 

0.046 

4.82 

ml8 

10.7 

2.2 

24.9 

0.68 

1.78 

0.047 

1.8 

0.057 

4.8 

mpll 8 

2.45 

0.54 

35. 

0.3 

0.72 

2.57 

1.6 

0.035 

5.0 

Kecks 

- 

- 

- 

- 

- 

- 

2.3 

0.6 

7.56 

ml9 

0.003 

0.3 

27.3 

0.74 

20.7 

0.15 

2.5 

0.45 

8. 

mpll 9 

0.03 

0.095 

34.4 

0.14 

0.64 

1.95 

2.1 

0.05 

7.6 

Keck20 

- 

- 

- 

- 

- 

- 

0.096 

0.27 

5.9 

m20 

0. 

0.03 

23.8 

0.07 

14.0 

3.72 

0.11 

0.15 

6.4 

mpl20 

420. 

190. 

38. 

85. 

22.5 

7.9 

0.08 

0.43 

6. 


galaxies from the radio survey of Humphrey et al (Table 3) do 
not show the high Vs component which is neither present in the 
galaxies selected from the Vernet et al survey (Table 5). Here, 
the clouds are characterised by the very low Vs and no, similar 
or lower than those calculated to reproduce the Finkelstein et al 
spectra (Tables 8 and 9). 

We report in Fig. 8 Hell 1640/CIV 1549 versus Lya /CIV 
1549 calculated by the models presented in GRIDl, considering 


gas in the conditions suitable to the NLR of AGN. The line ra¬ 
tios were calculated by GRIDl adopting two different geometri¬ 
cal thickness of the emitting clouds, D=10^^ cm and lO'^ cm. In 
GRIDl, for each set of no and Vs a series of models with differ¬ 
ent F, (log ^"=8,9,10,11,12) were adopted. F =0 represents the 
shock dominated case and corresponds to the lowest Hell/CIV 
and Lya /CIV. Solar (Allen 1976) element abundances were used 
to calculate the spectra. The Hell/CIV and Lya /CIV trends (for 


14 



M. Contini: Lyo- and UV line models for galaxies at z<3.7 



Fig. 8. The Hell/CIV versus Ly/CIV for NLR of AGN. Open circle : Norman et al, hexagon : Dawson et al, diamond ; Stern et al, 
large pentagram ; Dey et al; red asterisks ; Collins et al., triangles : Vernet et al, open squares ; Humphrey et al.;filled pentagrams 
(small for SB and large for AGN models):Ly et al. ; black filled squares (small for SB and large for AGN models) : Finkelstein et 
al.; small circle: Erb et al; plus: Fosbury et al. Line symbols : see Table 1. 


each Vs ) follow the increasing photoionization flux, because 
both Hell and Lya are recombination lines which depend on 
photoionization, while the fractional abundance distribution of 
carbon ions in different ionization levels strongly depends on 
gas temperature and density which are affected by the shock. 

The data observed for nearly all the AGN dominated objects 
are nested between models showing Vj ~ 100 and 200 km s ', 
the fluxes ranging throughout a large scale, and both geometri¬ 
cally thin and thick clouds. Also shown in Fig. 8 are the spectral 
data observed from the local galaxy Mrk3 (Collins et al 2005) 
because they include the Lya and UV lines. The data (reddening 
corrected) refer to different regions within the galaxy. The spec¬ 
tra represent average conditions covering small regions within 
the galaxy. So they tend to show intrinsic conditions of the gas 
rather than the averaged ones obtained by modelling the ob¬ 
served spectra covering the entire galaxy, as necessarily occurs 
for high z galaxies. The Mrk3 data are less scattered than those 
observed from high z galaxies and follow the model trends. 

In Fig. 8 we added the data based on the Finkelstein et 
al and Ly et al. observations. They do not represent observed 
Hell/CIV versus observed Lya /CIV, but the calculated line ra¬ 
tios which are strongly constrained by the observed optical lines. 
Finkelstein et al and Ly et al. objects follow the trends of AGN 
(H-shock) dominated models with roughly constant velocities. 
The observed line ratios referring to the other galaxies are scat¬ 


tered for different F. The Finkelstein et al and Ly et al samples 
are situated close to higher shock velocity models. 

Finally, in Figs. 9 the physical parameters (Vj , no, T, , 
U, F, D), the Ha flux at the nebula and O/H and N/H relative 
abundances calculated for the present sample of galaxies are 
presented as function of the SLR. SLR were calculated by the 
Ha flux measured at Earth by the observers. We chosed Ha in¬ 
stead of Lya because less affected by scattering and absorption 
problems. For galaxies in which the Ha data were not available, 
we used Lya and we calculated Ha line fluxes by noticing that 
the Lya /H/3 line ratios can be approximated to ~25 for many 
models (Fig. 1) and that Ha /H/? ~ 3. 

The sample of galaxies which emit Lya and UV lines do not 
show peculiar physical conditions of the emitting gas, nor ab¬ 
normal element abundances. However, for some galaxies of the 
Humphrey et al, Dawson et al. Stern et al, Dey et al, Norman et 
al surveys, emitting gas with relatively high Vs (> 1000 kms^*) 
is revealed by the complex FWHM of the line profiles. Such 
velocities are upper limits to those characteristic of the AGN 
NLR. Generally, higher no accompany higher Vs leading to col- 
lisionally deexcited forbidden lines. This is not the case for 
the Lya line emitting galaxies modelled throughout this paper. 
In the Humphrey et al galaxy 0211-122 spectrum reported in 
Table 2, for instance, [OIII] 5007 lines with FWHM = 1300+60 
kms-'and [Oil] 3727 with FWHM=730±160 kms-'were ob¬ 
served. Most of the high velocity models are shock dominated. 
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In these cases, the electron temperature and density steeply drop 
to T>10^ K (Fig. 6) at a certain distance from the shock-front 
(which depends on Vs ) because the primary flux is absent and 
the secondary one is relatively weak. 

The high T* (> 10^ K) which indicates that the SB stars are 
close to outburst (Contini 2014b) were found in a few objects 
(4C-1-23.56 at z=2.479 from the Humphrey et al survey, SST24 
at z=2.858 from the Dey et al survey and the Lynx arc at z=3.357 
from Fosbury et al observations ) within more than 50 galaxies 
in the sample. These outstanding temperature, which are charac¬ 
teristic of galaxies showing some activity (Contini 2014b), are 
isolated events at the state of the art. Considering the evolution¬ 
ary paths in the H-R diagrams for stars of different mass during 
the pre-main-sequence phase, these temperatures can be reached 
by >9Mo stars in a lifetime of > 10^ years (Prialnik 2000, fig. 
8.1), which is very short compared to the redshift z age (> 2 
Gyr). 

Fig. 9 shows that the SFR between ~ 1 and ~ 300 M© 
yr“'is accompanied by relatively high shock velocities (> 1000 
kms“^), but gaseous cloud densities are relatively low, similar 
to those in the AGN NLR and in the SB surrounding (no<350 
cm-3) for SFR > 10 Mq yr *. O/H reaches a lower limit min¬ 
imum for some galaxies at SFR< 10 M© yr“', then it grows 
with SFR to solar values at SFR ~ 300 M© yr“^ The N/H rel¬ 
ative abundances show a large scattering. N/H are particularly 
low for a few galaxies in the Finkelstein et al survey at z~ 0.25. 
Humphrey et al, however, claim that nuclear activity and star for¬ 
mation which lead to solar metallicities at early times, are seen 
through lack of variability and high metallicity between z=2.5 
and 0. The geometrical thickness of the emitting clouds presents 
minimum D < 10*® cm. D and consequently also Ha fluxes cal¬ 
culated at the nebula show a large scattering covering roughly 4 
orders of magnitude for SFR between 1 and 650 M© yr *, indi¬ 
cating fragmentation by turbulence in a shock dominated regime. 
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Fig. 9. The physical parameters calculated for the galaxy surveys as a function of the SFR. Symbols as in Contini 2014a, hg. 5. For 
the present sample ; blue open squares : Finkelstein et al for AGN; red open squares : Finkelstein et al. for SB. The results for Ly et 
al (2014) data are represented by pentagrams, red for starbursts and black for AGN; black hlled circles: Dey et al., Dawson et al., 
Norman et al.. Stern et al., Humphrey et al., Vernet et al., Erb et al, and Fosbury et al. AGN models; red hlled circles : the same for 
SB dominated models. 
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